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1.0  Scope. 

This  document  describes  the  construction,  deployment,  and  initial  data  processing 
of  four  optical  drifting  buoys.  The  buoys,  which  measure  spectral  water-leaving 
radiance  and  downward  irradiance,  were  successfully  deployed  in 
October/November,  1993  in  the  Equatorial  Pacific  in  collaboration  with  the  iron- 
seeding  experiment  IRONEX. 

2.0.  Optical  Drifting  Buoys. 

The  Ocean  Color  Monitor  (OCM)  optical  drifting  buoy  represents  a  novel  means 
to  autonomously  acquire,  and  retrieve,  data  on  upwelling  spectral  radiance  ("ocean 
color")  and  surface  downwelling  irradiance  from  remote  ocean  locations.  The 
buoys  are  designed  to  be  either  ship,  or  aircraft,  launched.  Data  is  automatically 
telemetered  in  near  to  real  time  via  CLS/ARGOS  and  Internet  to  Dalhousie 
computers.  The  buoys  also  measure  sea-surface  temperature. 

Further  details  of  the  buoy  construction,  optical  characterization,  calibration  and 
data  reporting  can  be  found  in  Appendix  A,  Ocean  Color  Monitor  User's  Manual. 

3.0.  Deployment  and  Operations 

The  buoys  were  packed  at  the  factory,  and  shipped  to  Miami  where  they  were 
loaded  on  the  Columbus  Iselin  to  take  to  the  deployment  site.  Instructions  were 
provided  to  the  ships  crew  on  deployment  techniques. 

The  buoys  were  deployed  in  association  with  the  joint  ONR/NASA/NSR  iron 
enrichment  study  south  of  the  Galapagos  Islands  in  October/November,  1993 
(IRONEX).  An  area  approximately  lOOkm^  was  fertilized  with  an  iron  sulphate 
solution  to  evaluate  the  effect  on  the  phytoplankton  population. 

After  the  patch  was  inoculated,  the  first  buoy  was  deployed  at  the  nominal  patch 
center  at  90oW  5°S  on  28  October,  1993.  The  next  two  subsequent  buoys  were 
deployed  two  days  and  6  days  later  respectively,  near  where  the  patch  center  was 
thought  to  be.  The  final  buoy  was  deployed  later  during  a  survey  of  the  Galapagos 
Island  plume  region  near  the  Equator. 

All  buoys  turned  on  and  began  transmitting  soon  after  deployment. 
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4.0  Data  Processing 

The  raw  ARGOS  telemetry  data  (Level  0  data;  file  extension  *.ARG)  is 
automatically  posted  once  per  day  to  Dalhousie  computers  via  Internet.  This 
hexidecimal  data  is  processed  to  a  Level  1  data  set  which  consists  of  earth-located, 
time-stamped  (time  of  transmission)  calibrated  radiances,  irradiance,  and 
temperature  (file  extension  *.DAT). 

The  Level  1  data  is  subjected  to  fiuther  processing  which  does  a  number  of 
operations.  First,  the  time  is  corrected  to  the  midpoint  of  the  hourly  averaged  data 
in  GMT.  Second,  the  calibrated  radiance  data  is  used  to  estimate  a  spectral 
attenuation  coefficient  to  propagate  the  radiance  to  and  through  the  sea-surface. 
Using  the  downwelling  irradiance  at  490  ran,  and  a  spectral  atmospheric 
transmission  model,  each  radiance  waveband  is  normalized  to  the  ratio  of  the 
spectral  irradiance  at  the  sea-surface  relative  to  the  extraterrestrial  atmospheric 
irradiance  to  produce  normalized,  water-leaving  radiances.  These  values  represent 
the  expected  water-leaving  radiances  with  the  sun  at  zenith  and  with  no 
intervening  atmosphere;  it  provides  a  means  to  compare  data  from  different 
illumination  conditions  as  well  as  a  direct  comparison  with  similar  quantities 
computed  from  satellite  (SeaWiFS)  observations. 

The  data  often  has  transmission  errors  (see  Users  Manual),  which  to  first  order  are 
corrected  at  this  stage.  Correction  methods  look  for  drops  in  battery  voltage,  and 
anomalous  temperature  readings.  Duplicate  transmissions  are  also  removed  at  this 
step.  Some  hand  processing  is  usually  required  to  remove  bad  data  points  not 
detected  by  the  processing  program. 

Finally,  a  first  order  pigment  concentration  is  computed  from  the  normalized 
water-leaving  radiances  using  the  global  CZCS  algorithm.  For  low  irradiance 
levels  (i.e.  night-time),  the  last  good  pigment  estimation  is  carried  until  sufficient 
irradiance  is  available  to  provide  good  radiance  measurements. 

The  result  of  these  operations  is  a  Level  2  data  set  (file  extension  ♦.PRO).  These 
data  can  be  further  processed  if  desired  into  a  variety  of  averaged  or  experimental 
products. 

5.0.  Initial  Results  and  Discussion 

As  of  the  28th  of  December,  1993,  all  buoys  continued  to  transmit  data.  Buoys 
drifted  initially  in  a  northeasterly  direction,  then  almost  due  north.  Buoys  4742 
and  4743  were  deployed  very  close  to  each  other,  and  remained  close  throughout 
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the  past  two  months.  As  of  28  December,  they  were  within  30  km  of  each  other, 
despite  the  2000  km  they  have  translated  since  deployment.  All  buoys  showed  a 
high  degree  of  coherence  in  their  relative  velocities  and  direction  of  drift.  Figure  1 
shows  the  buoy  tracks  until  22  December,  1993. 

Figure  1.  Optical  Drifting  Buoy  Tracks. 

Buoy  Position  Tracks,  IRONEX 


Longitude 
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5. 1.  Sea-surrface  temperature.  The  data  on  sea-surface  temperature  shows  both 
the  long  term  variations  as  the  buoys  transited  over  their  range,  and  also  show  a 
clear  diurnal  SST  variation.  The  SST  plots  for  the  four  buoys  on  a  consistent  time 
axis  are  shown  in  Figures  2-5  below. 
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Figure  2.  Sea-sutface  temperature  variations.  Buoy  4745. 
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e  5.  Sea-surface  temperature  variations.  Buoy  4744. 
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Figure  4.  Sea-surface  temperature  variations.  Buoy  4743. 
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Figure  5.  Sea-surface  temperature  variations.  Buoy  4742. 
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5.2  Surface  Irradiance.  Surface  irradiance  at  490  imi  is  shown  in  the  Figures  6-9. 


Figure  6.  Surface  irradiance  (490  nm)  variations.  Buoy  4745. 
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Figure  7.  Surface  irradiance  (490  nm)  variations.  Buoy  4744. 
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5.3.  Surface  Pigment.  Surface  pigmr'it,  computed  with  the  nominal  CZCS 
algorithms  is  given  in  Figures  K  13  below. 

Figure  10.  Surface  pigment  variations.  Buoy  4745 
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Figure  12.  Surface  pigment  variations.  Buoy  4743 
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5.4.  Evaluation..  The  buoy  performance  was  nominal  for  the  duration  to  date. 
SST's  and  pigment  estimates  are  comparable  to  those  measured  during  the 
IRONEX  ship  observations  and  are  consistent  with  local  climatology.  Of 
particular  note  are  the  pigment  time-series  for  buoys  4742-4745  which  show  a 
high  degree  of  correspondence  in  time.  This  is  the  first  time  buoys  have  been 
deployed  which  can  be  considered  essentially  replicaties  and  these  results  are 
encouraging. 

6.0  Data  Availability. 

Level  2  data  is  available  for  anonymous  ftp  on  /pub/ironex  on 
predator.ocean.dal.ca.  Those  wishing  raw  or  calibrated  radiance  should  contact 
marlon@predator.ocean.dal.ca  and  it  can  be  made  available  as  well. 
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INTRODUCTION 

The  OCM  (Ocean  Color  Motutor)  is  designed  to  be  a  self-contained  precision 
optical  instrument  suitable  for  expendable  launch  in  harsh  ocean  tenons  to  remotely 
monitor  the  apparent  optical  properties  of  the  upper  ocean.  The  system  can  be  launched 
from  multiple  platforms,  including  aircraft  and  ships  of  opportunity,  and  relays  data  back 
to  the  user  via  the  ARGOS  satellite  system. 

The  OCM  has  a  seven-channel  upwelling  radiance  sensor  and  a  single  channel 
downwelling  irradiance  sensor.  The  system  uses  a  proprietary  filter/photodiode  system  to 
provide  improved  signal  performance,  ruggedness  and  sensor  stability  necessary  for 
oceanographic  use.  The  OCM,  designed  as  a  multiplatform  module,  has  been  used 
successfully  as  a  ship-launched  drifting  buoy,  an  air-launched  drifting  buoy,  moored 
systems,  and  tethered  systems. 

The  rugged  design  of  the  OCM  buoy  resulted  from  its  development  as  an  aircraft 
gravity  tube  launch  system  which  was  required  to  withstand  impacts  on  the  ocean  surface 
at  spe^s  of  over  50m/sec.  This  was  verified  through  test  launches  from  NASA's  P3B 
aircraft.  The  OCM  has  also  been  field  tested  in  all  deployment  modes  in  temperatures 
ranging  from  -20®C  to  +40®C  with  no  degradation  in  performance. 

The  filters  used  in  the  OCM  were  custom  designed  to  match  the  full  20nm 
bandwidth  of  the  upcoming  SeaWiFS  satellite  with  an  additional  lOnm  channel  at  683nm 
for  the  measurement  of  solar-stimulated  fluorescence.  The  use  of  custom  filters  allows  us 
to  precision  match  instruments  for  our  customers.  The  spectral  scans  for  these  filters  are 
available  upon  request. 

A  significant  amount  of  design  and  field  testing  has  gone  into  the  OCM  since  the 
first  prototype  in  1989  and  we  are  confident  you  will  be  pleased  with  its  performance  in 
any  ocean  environment. 


Scott  McLean, 
Project  Engineer 


Telephone:  (902}  492-1780 


Facsimile:  (902)  492-4781 


Satlantic  Inc. 


Richmond  Terminal,  Pier  9 
3295  Barrington  Street 
Halifax,  Nova  Scotia 
Canada  B3K  5X8 


1.0  General  Specifications 


Upwelling  Radiance  Sensor  Characteristics 


Sensor  Model:  OCR- 100 


Spatial  Characteristics: 


-  Field  of  view:  5®  (0.025  steradians)  in  water 

T  (0.100  steradians)  in  air 

-  Entrance  aperture;  4.78  mm  diameter 

-  Detectors;  custom  13  mm^  silicon  photodiodes 


Spectral  Characteristics: 


-  Bandwidth  range;  400-700  nm 

-  Number  of  channels;  7 

-  Spectral  bandwidth:  6  channels  20  nm 

1  channel  10  nm 

-  Filter  Type,  custom  low  fluorescence  interference 

-  Discrete  wavelengths  (centers):  412,  443,  490,  510,  555,  670,  683  nm 


Optical  Characteristics: 


Out  of  band  rejection;  10"^ 
Out  of  field  rejection;  5x10’^ 


Temporal  Characteristics: 


-  System  time  constant;  0.015  seconds 

-  -3dB  frequency:  10  Hz 
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Downwelling  Irradiance  Sensor  Characteristics 


Sensor  Model :  ED-100 


Spatial  Characteristics: 


-  Field  of  view:  cosine  response 

-  Collector  area:  86.0  mm^ 

-  Detectors:  custom  13  mm^  silicon  photodiodes 


Spectral  Characteristics: 


-  Bandwidth  range:  400-700  nm 

-  Number  of  channels:  1 

-  Spectral  bandwidth:  20  nm 

-  Filter  Type:  custom  low  fluorescence  interference 

-  Discrete  wavelengths  (center):  490nm 


Optical  Characteristics: 


-  Out  of  band  rejection:  10'^ 

-  Cosine  response:  within  3%  3-60° 

within  10%  60-89° 


Temporal  Characteristics: 


-  System  time  constant:  0.015  seconds 

-  -3dB  frequency:  10  Hz 


Operational  Characteristics:  (Buoy  System) 

-  size:  91cm  X  12.Scm  diameter  (packaged  air-launch  only) 
-weight:  12kg 

-  power:  internal  -  50  alkaline  'C  cells 

-  operating  lifetime:  up  to  6  months 

-  temperature  rating:  -20°C  to  +40°C  (operational) 

-  satellite  transmitter:  MetOcean  MAT8602  ARGOS  PTT 
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2.0  Absolute  Radiometric  Calibration  (Radiance) 

Absolute  radiometic  radiance  calibration  is  done  using  a  calibrated  lOOOW  FEL 
lamp  on  a  Sm  optical  bar  using  the  'plaque  method*.  The  lamp  is  powered  by  an  Optronics 
83DS  current  source.  The  flux  from  the  lamp  is  normally  incident  on  a  50  cm  diffuse 
reflectance  target  standard  at  a  distance  of 200.0  cm.  The  instrument  views  the  target  at 
an  angle  of  45. 0  deg  such  that  the  field  of  view  of  all  the  sensors  is  completely  covered  by 
the  target.  The  calibration  radiances  are  determined  using  equation  1  below; 

(1 )  L(;t)  =  (E(X,  50cm)  /  n)  *  (50.0  cm  /  200.0  cm)2  *  p(X) 

where: 

L(>.)  is  the  calibration  radiance 
E(7.,  SOctn)  is  the  lamp  calibration  at  SOcm 
(50.0  cni/200.0  cm)^  is  the  1/R^  distance 
P().)  is  the  reflectance  target  calibration 

Reflection  Target:  Labsphere  SRT-99-180  S/N  001873 
Standard  Lamp;  Hoffman  S/N  91615 
Voltmeter;  HP34401A  S/N  3146 A09840 


The  voltage  output  of  the  instrument  is  measured  on  an  HP34401 A  6.5  digit 
multimeter  at  a  sample  rate  of  20Hz.  A  5  second  average  is  used  for  the  calibration 
output. 


wavelength 

(nm) 

lamp  iiradiance 
(^W/cm^  /nm) 

target  reflectance  calibration  radiance 
(pW/cm^  /nm/sr) 

OCM  input 
(mV) 

OCM  dark 
(HV) 

410 

2.766 

0.981 

0.0540 

86.45 

-f610 

444 

4.388 

0.982 

0.0857 

115.12 

+280 

489 

7.219 

0.981 

0.1409 

283.10 

-10 

511 

8.671 

0.982 

0.1694 

298.80 

-60 

553 

11.55 

0.985 

0.2263 

476.70 

+740 

668 

18.60 

0.982 

0.3634 

1372.40 

+660 

684 

19.36 

0.983 

0.3786 

1274.80 

•640 

Table  I  -  OCM  ID  4742  Absolute  Radiometric  Calibration  (Radiance) 
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2.1  Absolute  Radiometric  Calibration  (Irradiance) 

Absolute  radiometic  irradiance  calibration  is  done  using  a  calibrated  lOOOW  FEL 
or  DXW  lamp  on  a  5m  optical  bar  using  direct  radiation  from  the  lamp.  The  lamp  is 
powered  by  an  Optronics  83DS  current  source.  The  flux  from  the  lamp  is  normally 
incident  on  the  irradiance  sensor  cosine  collector  at  a  distance  of  100.0  cm.  The 
calibration  irradiances  are  determined  using  equation  2  below. 

(2)  E(X,  75cm)  =  E(X,  50cm)  •  (50.0  cm  /  75.0  cm)2 
where: 

£(>..  7Scin)  is  the  calibration  irradiance 
£().,  SOctn)  is  the  lamp  calibration  at  50cm 
(50.0  cm/75.0  cm)^  is  the  l/R^  distance 

Standard  Lamp.  Hoffman  S/N  91604 
Voltmeter:  HP34401A  S/N  3146A09840 

The  voltage  output  of  the  instrument  is  measured  on  an  HP34401 A  6.5  digit 
multimeter  at  a  sample  rate  of  20Hz.  A  5  second  average  is  used  for  the  calibration 
output. 


wavelength 

lamp  irradiance 

calibration  irradiance  OCM  input 

OCM  dark 

(nm) 

(pW/cm^/nm) 

(uW/cm^/nm)  (mV) 

(UV) 

489 

7.219 

3.208  50.% 

+10 

Table  2  -  OCM  ID  4742  Absolute  Radiometric  Calibration  (Irradiance) 
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2.2  Immersion  Effect  (Radiance) 

Due  to  the  difference  in  indices  of  refraction  between  air  (where  the  instrument  is 
calibrated)  and  water  (where  it  is  operated)  a  correction  factor  must  be  applied  to  obtain 
the  effective  in  water  radiances.  This  correction  factor  is  referred  to  as  the  immersion 
factor.  There  are  two  effects  contributing.  First,  the  reduction  in  solid  angle  viewed  by 
the  sensors  effectively  reduces  the  amount  of  flux  into  the  sensor.  This  correction  is  given 
by  FI. 

(2)  FI  (X.)  =  (ti^X))2  where  n*  die  inde.x  of  refraction  of  water 


To  correct  for  calibration  values  in  air,  the  in-water  values  are  multiplied  by  the  effective 
loss  of  viewing  area  in  water  (FI ). 

The  second  effect  is  due  to  the  change  in  index  of  refraction  at  the  glass/air 
(glass/water)  interface.  This  correction  is  given  by  F2. 


(3)  F2(X)  =  (nv^X)  +  ilg(>.))2  /  ((T|*(X)  •  (I  + 

where  T|g  is  the  index  of  refraction  of  window 


Since  the  indices  of  refraction  of  water  and  glass  are  better  matched,  there  are  less 
reflection  losses  at  the  window.  The  immersion  factor  thus  reduces  the  in-water  values  to 
correct  for  this  effect. 

The  total  immersion  effect  is  then; 

(4)  Imm(X)  =  F2(X)*F2(X) 

Thus  the  correction  for  actual  in-water  radiance  values  is; 

(5)  Lwater  actuai(X)  =»  Lwater  measured(X)  •  Imm(X) 
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The  results  for  the  OCM  are  in  table  3  below: 


wavelength 

index  water 

index  window 

immersion 

(nm) 

correction 

410 

1.349 

1.534 

1.747 

444 

1.347 

1.534 

1.741 

489 

1.344 

1.529 

1.734 

511 

1.343 

1.529 

1.732 

553 

1.341 

1.525 

1.728 

668 

1.338 

1.520 

1.720 

684 

1.337 

1.520 

1.719 

Table  3  -  OCM  Immersion  Corrections 


References: 

Austin,  R.W.  (1976),  Air-water  radiance  calibration  factor.  Tech.  Memo. 

ML-76-004t,  Vis.  Lab.,  Scripps  Institute  of  Oceanography,  8pp. 
Austin,  R.W.,  and  G.  Halikas  (1976),  The  index  of  refraction  of  seawater, 
SIO  Ref  76-1,  Vis.  Lab.,  Scripps  Institute  of  Oceanography,  64pp. 
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3.0  Calibration  Summary  OCM  ID  4742 

The  CMTM  buoy  transmits  the  data  from  the  sensors  in  hexadecimal  counts  from  00 
to  FF  (0  to  2SS).  The  minimum  voltage  that  can  be  read  is  1  SmV  (one  count),  the 
maximum  voltage  that  can  be  read  is  3.82SV  (25S  counts).  The  calibration  summary  will 
give  the  results  in  counts  such  that  the  calibration  factor  (F)  can  be  applied  as  in  equations 
(1)  and  (2)  below; 

( 1 )  radiance  =  counts  •  F  «  immersion  (fiW/cm^/nm/sr) 

where  immersion  comes  from  Table  3 

(2)  irradiance  =  counts  •  F  (jiW/cm^/nm) 

The  calibration  factor  (F)  is  computed  using  equation  (3) 


(3)  F  (units/count)  =  calibration  (units)  /  calibration  volts  (mV)  •  1 S  (mV/count) 


Sensor 

Filter# 

Lot# 

CWL 

(nm) 

Calibration 

(mV) 

Cal  Radiance 
(jiW/cm^/nm/sr) 

F 

Saturation* 

(^W/cm^/nm/sr) 

Lu412 

25 

1493 

410 

85.84 

0.0540 

9.436E-3 

2.406 

Lu443 

24 

1293 

444 

114.84 

0.0857 

11.19E-3 

2.854 

Lu490 

17 

1093 

489 

283.11 

0.1409 

7.465E-3 

1.904 

Lu510 

19 

1393 

511 

298.86 

0.1694 

8.502E-3 

2.168 

Lu555 

29 

1393 

553 

475.96 

0.2263 

7.132E-3 

1.819 

Lu670 

16 

1393 

668 

1371.7 

0.3634 

3.974E-3 

1.013 

Lu683 

23 

1393 

684 

1275.4 

0.3786 

4.453E-3 

1.135 

^saturation  is  calculated  at  2S5  counts 

Table  4 

OCR-1 00  radiance  sensor  S/N  055 

Sensor 

Filter# 

Lot# 

CWL 

(nm) 

Calibration 

(mV) 

Cal  Irradiance 
(nW/cm^/nm) 

F 

Saturation* 

(jiW/cm^/nm) 

Ed490 

18 

1093 

489 

50.95 

3.208 

0.944 

240.8 

*saturation  is  calculated  at  235  counts 

Table  5  ED-100  irradiance  sensor  S/N  055 
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This  is  a  copy  of  the  calibration  file  04742A.CAL  for  use  with  the  0CM_PR03 
and  0CM_PR04  processing  programs.  It  contains  the  F  values  used  in  equation  (3) 
above  to  calculate  physical  units  from  the  buoy  digital  counts. 

#  Calibration  File  -  OCM  Ship-Launch  ID4742 
#OCR-100S/N055  ED-100  S/N  055 

#  Calibration  Date;  19  September  1993 

#  order;  Lu683,  Lu670,  Lu555.  Lu510,  Lu490.  Lu443,  Lu412.  Ed490 

# 

4.453E-3  0 
3.974E-3  0 
7.132E.3  0 
8.502E-3  0 
7.465E-3  0 
11.19E-3  0 
9.436E.3  0 
0.944  0 
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4.0  Deployment 

The  OCM  has  three  standard  launch  modes;  air-launch  deployment,  ship-launch 
airbag  deployment,  and  solid  collar  deployment.  Air-launch  and  ship-launch  airbag 
deployments  are  by  far  the  simplest  as  the  package  is  designed  to  deploy  autonuitically. 
These  packages  have  salt  water  activated  squibs  which  automatically  deploy  the  buoy 
when  in  the  water.  The  airbag  floatation  collar  on  these  packages,  however,  is  not  as 
robust  as  the  solid  collar.  The  solid  collar  units  are  designed  for  experiments  in  extremely 
harsh  conditions  (such  as  ice)  or  in  moored  configurations  where  the  float  may  chaff  on 
mooring  lines. 

OCM  ARGOS  ID  #4742  is  a  Ship-Launch  Solid  Collar  Buoy.  See  Section  4.3 
for  deployment  details. 


4.1  Air-Launch  Package  Deployment 

The  OCM  has  a  US  Navy  NAVOCEANO  air-launch  certification  as  a  CMOD  type 
buoy  and  is  suitable  for  gravity  tube  air  launch  as  an  A-size  sonobuoy  package.  It  has 
been  sucessililly  air-launched  from  NASA's  P3B  aircraft  with  a  100%  sucess  rate. 
Although  the  package  is  certified  over  a  large  window  of  deployment  speeds  and  altitudes, 
a  speed  of  200  knots  (or  less)  and  an  altitude  of  at  least  2000  feet  is  recommended  to 
reduce  the  probability  of  package  damage  on  impact  with  the  water.  This  altitude  assures 
that  the  horizontal  velocity  of  the  buoy  is  near  zero  at  impact  and  is  in  a  near  vertical 
position. 

The  standard  air-launch  package  is  shipped  in  an  SLC  (sonobuoy  launch  container 
-  13cm  diameter  by  100cm  long)  which  looks  much  like  an  AXBT  package  (but  heavier). 
The  ARGOS  ID  for  the  OCM  inside  the  conatiner  is  written  on  the  lid  of  the  SLC  as  well 
as  on  the  side.  Do  not  store  the  OCM  outside  of  the  SLC  before  launch,  as  damage  to  the 
parachute  may  occur. 

To  launch  the  buoy,  remove  the  retaining  clip  holding  the  lid  on  the  SLC  and  pry 
off  the  SLC  lid.  Slide  the  buoy  out.  The  ARGOS  ID  is  also  written  on  top  of  the  buoy's 
wind  flap  which  opens  its  parachute  on  launch.  When  the  pilot  indicates  it  is  safe  to 
launch,  insert  the  buoy  into  the  launch  tube  wind  flap  end  first  (failure  to  do  this  may 
damage  the  aircraft!)  and  drop  the  buoy  out  of  the  launch  tube. 
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4.2  Ship-Launch  Airbag  Package  Deployment 

The  ship-launch  airbag  OCM  is  designed  to  be  an  easy  self-deploying  package  that 
can  be  deployed  from  any  vessel  moving  at  full  speed.  It  is  very  similar  to  the  air-launch 
package  except  that  it  does  not  have  a  parachute  and  windflap. 

The  standard  ship-launch  package  is  shipped  in  an  SLC  (sonobuoy  launch 
container  -  13cm  diameter  by  100cm  long).  The  ARGOS  ID  for  the  OCM  inside  the 
container  is  written  on  the  lid  of  the  SLC  as  well  as  on  the  side.  The  side  of  the  SLC  will 
also  be  clearly  marked  ‘ship-launch'  indicating  that  the  unit  does  not  have  a  parachute. 

To  launch  the  buoy,  remove  the  retaining  clip  holding  the  lid  on  the  SLC  and  pry 
off  the  SLC  lid.  Slide  the  buoy  out  The  buoy  cannot  be  deployed  inside  the  SLC. 

The  ARGOS  ID  is  also  written  on  top  of  the  clam  shell  retainer.  The  buoy  can  be  tossed 
over  the  side  of  the  vessel  for  deployment.  Do  not  get  the  buoy  wet  on  deck  as  it  may 
activate  the  squib,  blowing  the  clam  shells  off  and  injuring  personnel.  Note  that  the 
buoy  will  descend  to  a  depth  of  10m  before  the  squib  fires,  which  takes  about  10-15 
seconds.  The  squib  will  fire,  inflating  the  airbag,  blowing  off  the  clam  shells,  releasing  the 
drouge  and  deploying  the  antenna.  The  buoy  will  then  rise  to  the  surface  and  begin 
transmitting.  Do  not  launch  the  buoy  from  a  stationary  vessel  as  the  buoy  may  strike 
the  vessel  on  its  way  back  to  the  surface,  damaging  the  antenna  and/or  the  cosine 
collector. 

4.3  Ship-Launch  Solid  Collar  Deployment 

This  OCM  has  a  high  density  foam  solid  collar  which  can  withstand  very  harsh 
conditions  without  being  damaged.  Since  the  collar  is  a  fixed  size,  the  package  is  not  as 
easily  deployed  as  the  airbag  type,  but  is  more  suitable  for  moored  type  applications. 

The  buoy  comes  packaged  in  a  37cm  diameter  tube,  100cm  long.  The  package 
does  not  automatically  deploy  and  must  be  deployed  manually.  The  top  cover  of  the  tube 
has  the  buoy  ARGOS  ID  written  on  it.  Remove  the  screws  which  secure  both  covers  and 
slide  the  packing  tube  off  the  buoy.  The  ARGOS  ID  is  on  a  stamped  plate  which  is 
attached  to  the  outer  casing  (drouge),  just  above  the  magnet.  The  magnet,  when 
removed,  turns  on  the  ARGOS  transmitter.  Remove  the  restraint  holding  down  the 
antenna,  exposing  the  irradiance  cosine  collector.  Pull  the  antenna  to  its  fullest  extent. 
Remove  tape  holding  the  drouge  (outer  casing)  on  and  slide  it  off  (this  also  turns  the 
transmitter  on  as  the  magnet  slides  off).  The  drouge  can  is  attached  to  the  buoy  via  a  20m 
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long  cable  to  improve  its  stability.  The  drouge  cable  is  packaged  onto  the  side  of  the 
radiance  sensor  with  wax  to  prevent  tangling.  With  the  vessel  stopped,  lower  the  drouge 

into  the  water  (if  the  vessel  is  not  stopped  the  drouge  acts  as  a  sea  anchor  and  will  pull  the 
buoy  off  the  deck,  into  the  water,  possibly  damaging  it).  When  the  drouge  is  fully 
deployed  the  buoy  can  be  lowered  or  dropped  vertically  into  the  water.  Care  must  be 
taken  to  not  break  the  antenna. 

Mooring  schemes  can  be  used  effectively  with  this  buoy  by  tightly  clamping  a  band 
around  the  center  of  the  float  which  contains  attachment  points.  The  band  should  dig  into 
the  float  such  that  it  cannot  slide  off.  The  most  effective  mooring  for  estuary  work  is  a 
three  point  mooring  which,  with  three  guard  floats  arranged  in  a  triangle  about  10m  on  a 
side.  The  buoy  can  then  be  tied  off  to  each  float,  providing  a  very  stable  system.  The 
buoys  have  also  been  moored  to  a  single  float,  but  this  destabilizes  the  buoy  causing 
amplified  wave  induced  tilts,  increasing  data  loss  to  to  poor  transmission  (see  Data 
Processing  section  for  details).  Many  other  mooring  schemes  are  possible,  contact  us  to 
discuss  your  application. 
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5.0  Data  Processing 

ARGOS  data  is  acquired  from  the  ARGOS  Data  Processing  Center  in  Landover, 
USA  via  the  TYMNET  system.  The  raw  ARGOS  is  downloaded  in  the  TX  format  (see 
ARGOS  Users  Manual  section  3.2.5)  to  be  compatible  with  the  data  processing  software. 
The  data  consists  of  sets  of  data  points  referred  to  as  'hits'.  Each  'hit'  is  a  data  point  sent  by 
the  buoy  that  was  stored  by  a  satellite-based  ARGOS  receiver.  The  buoy  transmits  data  at 
90  second  intervals  continuously.  Reception  at  the  satellite  is  determined  by  the  buoy's 
position  and  the  satellite  orbit  (typically  10-12  times  per  day).  The  data  is  then  sent  by  the 
satellite  to  the  next  earth  station  it  passes  over,  where  it  is  then  processed  by  Service 
ARGOS  to  determine  the  platform's  position.  Raw  telemetry  (Level  0)  data  is 
downloaded  every  second  day  by  the  user  and  saved  with  an  .ARG  extension  A  sample 
data  'hit'  is  shown  below 


07188  44.672N  63 . 603W  0  343/1652Z-343/1651 

{  1)  2000468  301066F  301056C  301056E 

100025D  301046C  2000463  lOAO 

45  B6 

The  first  line  contains  the  platform  ID  (in  this  case  07188),  the  platform's  estimated 
position  in  decimal  degrees  (in  this  case  our  test  mooring),  the  location  accuracy  (single 
digit  from  0  to  3  -  0  means  accuracy  not  determinable,  I  within  1000m,  2  within  300m,  3 
within  150m)  and  then  the  day  and  time  of  the  data  collection  (in  Julian  day  and  time  in 
UTC).  The  second  day  and  time  is  the  last  position  update. 

The  number  in  brackets  (  1 )  on  the  second  line  is  a  compression  index  and  is 

ignored  by  data  the  processing. 

The  remaining  lines  contain  the  raw  telemetry  from  the  buoy  coded  in  hexadecimal. 
The  first  seven  groups  of  numbers  are  the  radiance  channels,  the  eighth  is  the  mast 
irradiance  sensor,  the  ninth  is  the  internal  thermistor  and  the  last  contains  the  battery 
voltage  and  the  tick  timer.  The  radiance  channel  number  contains  eight  hexadecimal  digits. 
If  not  all  digits  are  shown,  the  leading  blanks  are  zeros.  The  eight  di^ts  are  in  groups  of 
two  hexadecimal  digits  bytes).  The  first  two  digits  represents  the  channel  average.  This 
average  is  computed  by  taking  the  last  hour  of  data  (40  data  points)  and  computing  the 
average.  The  average  is  not  a  running  average;  the  previous  hour's  average  of  data  is 
transmitted  while  a  new  one  is  being  collected.  Every  hour  the  data  transmitted  is 
updated,  and  the  old  average  is  discarded.  The  second  two  digits  contain  the  minimum 
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value  recorded  during  the  hour.  The  third  two  digits  contain  the  maximum  value  recorded 
during  the  hour.  The  last  two  digits  contain  an  encoded  sum  of  squares  such  that  the 
standard  deviation  can  be  computed. 

In  order  of  occurrence,  the  radiance  channels  are  683nm.  670nm,  555nm,  5  lOnm, 
490nm,  443nm  and  412nm.  The  average,  minimum  and  maximum  are  converted  to 
radiances  using  the  equation; 

Radiance  =  F  •  (counts  -  offset)  •  immersion  (|.iW/cm-/nnt/sr) 

The  standard  deviation  is  computed  by  first  recovering  the  sum  of  squares  using  the 
equation; 


sum_squares  =  /  36.4)yjoooO 

The  standard  deviation  (in  volts)  is  then; 

Std_Dev_volts  *=  sqrt((sum_squares  -  n  *  average-)  /  (n-1)) 


In  radiance  units; 


Std_Dev_rad  =  Std_Dev_volts  *  255  /  3.825  *  F  ♦  immersion 

he  calibration  factors  (F)  are  stored  in  the  platform's  calibration  file  which  is  made 
up  of  the  platform  ID  with  a  .CAL  extension. 

The  irradiance  channel  is  the  last  group  of  digits  on  the  second  line  and  contains 
four  hexadecimal  digits  which  represent  the  irradiance  average  and  sum  of  squares  (two 
digits  each).  The  irradiance  sensor  on  all  platforms  is  an  Ed490  sensor .  The  sensor  values 
are  converted  to  physical  units  using  the  equation; 


Irradiance  =  F  *  (counts  -  offset)  (jiW/cm^/nm) 


The  first  number  on  the  last  line  is  the  buoy  hull  temperature.  This  is  converted 
from  counts  to  degrees  Celsius  using  the  following  equation; 

Temperature  ®  counts  ♦0.16-5 
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The  last  number  contains  the  battery  voltage  and  the  tick  timer. 

The  tick  timer,  stored  in  the  lower  6  bits,  indicates  how  many  minutes  it  has  been 
since  the  data  in  the  transmit  buffer  was  last  updated.  This  value  will  change  every 
transmission.  To  get  the  actual  time  for  the  middle  of  the  transmission  average,  take  the 
transmission  time,  subtract  the  tick  time  minutes,  and  then  subtract  30  minutes  more  to  get 
the  center  of  the  hour  period. 

The  upper  two  bits  indicate  the  battery  voltage.  Each  bit  represents  1.5  volts  above 
10.5  volts,  at  which  point  the  platform  transmitter  will  fail.  This  value  must  be  carefully 
examined  each  day  data  is  collected  if  the  buoy  is  to  be  recovered.  Once  the  batteries  fail, 
there  will  be  no  more  position  data.  For  quick  look  purposes  the  battery  voltage  can  be 
determined  by  looking  at  the  first  digit  of  the  last  number,  as  follows. 


hexadecimal 

binary 

voltage 

Fx 

llllxxxx 

15.0V 

Ex 

1 1 1 Oxxxx 

15.0V 

Dx 

1 lOlxxxx 

15.0V 

Cx 

1 lOOxxxx 

15.0V 

Bx 

101 Ixxxx 

13.5V 

Ax 

101 Oxxxx 

13.5V 

9x 

100 Ixxxx 

13.5V 

8x 

lOOOxxxx 

13.5V 

lx 

Oil  l.KXXX 

12.0V 

6x 

01  '-CxaXX 

12.0V 

5x 

OlOixxxx 

12.0  V 

4x 

OlOOxxxx 

12.0V 

3x,2x,lx,0x 

OOxxxxxx 

10.5V 

TX  Format  Battery  Voltage  Quick  Look  Table 
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5.0.1  Data  Transmission  Errors 

Due  to  the  256  bit  per  transmission  limit  on  ARGOS  there  is  no  checksum 
transmitted  with  the  data.  Erroneous  data  can  therefore  occur  at  random  and  must  be 
manually  removed.  Data  errors  will  occur  more  frequently  if  the  buoy  is  moored  or 
restrained  such  that  wave  action  causes  the  antenna  to  tilt  quickly  since  the  mooring  or 
restraint  will  not  follow  the  motion  of  the  buoy.  Best  performance  occurs  when  the  buoy 
is  allowed  to  free  drift  as  its  dynamics  are  designed  for  this  mode  of  operation. 

Examples  of  these  errors  (indicated  by  circled  data  records)  can  be  seen  in 
Appendix  3  where  0CM_PR04  has  been  run  on  file  SEP05.  ARG.  Buoy  04742  was 
moored  with  a  single  point  mooring  in  Monterey  Bay  in  50m  of  water  during  a  period  of 
strong  onshore  winds.  The  platform  antenna  was  moving  rapidly  during  this  period, 
resulting  in  poor  communications.  The  data  from  all  of  the  good  transmissions  shows  very 
consistent  data.  Buoy  04745  was  air-launched  into  the  equatorial  Pacific  and  was  free 
drifting;  it  occasionally  gets  a  bad  transmission. 

To  screen  bad  transmissions,  first  check  the  temperature.  Bad  transmissions  often 
have  very  high  or  very  low  temperatures.  Check  for  other  inconsistencies  such  as  very 
high  radiances  or  irradiances.  Bad  transmissions  will  be  obvious  when  you  plot  your  data. 
If  you  have  any  questions  about  this,  call  us. 
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S.l  Data  Processing  with  0CM_PR03.C 

Two  sample  programs  are  included  in  the  Appendix  1  0CM_PR03  .C  and 
0CM_PR04.C.  These  programs  read  TX  formatted  ARGOS  data  and  produce  calibrated 
ASCII  output  files  for  each  platform  transmitter.  These  programs  are  the  ones  that  we  use 
to  process  our  own  data,  and  are  provided  free  of  charge.  Please  feel  free  to  modify  them 
as  you  require.  We  use  0CM_PR03  .C  to  apply  calibrations  to  our  data  and  then  use 
MATLAB  to  produce  final  output  and  run  models. 

Here  is  a  brief  decription  of  0CM_PR03.C  Source  code  is  listed  in  Appendix  1 
and  examples  are  provided  in  Appendix  3: 


Program  0CM_PR03.C 

written  by;  Scott  McLean 
ver  date;  Jan  15/93 
distribution:  unlimited 
(c)  Satlantic  Inc.  1993 

valid  for  buoy  data  from  05  Aug  92  to  date 


This  program  ingests  OCM  type  buoy  datafiles  downloaded  from  the  ARGOS 
Processing  Center.  The  data  must  be  downloaded  from  the  processing  center  in  TX  format 
using  the  PRV  command  (refer  to  ARGOS  Users  Manual  section  3.2.5). 

The  program  extracts  the  buoy  data  from  these  files  and  produces  an  ASCII  file 
with  calibrated  data  in  physical  units.  Any  number  of  such  files  can  be  read.  The  data  in 
each  file  is  scanned  and  each  data  transmission  will  be  processed,  all  extraneous  characters 
(for  example,  the  login)  are  ignored.  Only  one  platform  can  be  processed  at  a  time,  if  you 
have  three  platforms,  you  must  run  the  program  three  times. 

This  program  produces  output  files  which  are  very  wide  (exceeding  256 
characters),  if  this  is  too  large  for  your  other  processing  programs,  use  0CM_PR04.C 
which  does  not  include  min,  max  and  standard  deviation  values  (only  the  averages)  and 
produces  a  more  compact  output  file. 


Telephone:  (902)  492-4780 


Facsimile:  (902)  492-4781 


Satlantic  Inc 


Richnwnd  Terminat,  Pier  9 
3295  Barrington  Street 
Halifax,  Nova  Scotia 
Canada  B3K  5X8 


Required  inputs: 

1)  INFILES.LST 

This  file  contains  a  list  of  all  files  that  are  to  be  processed.  The  program  assumes 
all  raw  data  files  have  the  .ARG  extension. 

2)  .ARG  files 

These  files  contain  the  data  obtained  on  various  days  from  the  ARGOS  processing 
center  in  TX  format.  All  files  must  have  the  .ARG  extension. 

3)  CAL  files 

These  files  contain  the  calibration  data  for  each  platform.  The  first  5  characters  of 
this  file  name  are  the  platform  ARGOS  ID  number.  The  program  will  prompt  you  for  the 
calibration  file  name.  It  will  remove  the  first  five  characters  from  the  calibration  file  name 
and  then  search  all  data  files  in  INFILES.LST  for  data  hits  containing  this  ID. 

4)  Output  file 

This  is  an  ASCII  output  file  containing  the  data  for  the  platform  requested  The 
program  will  prompt  you  for  a  file  name.  If  the  file  already  exists,  it  will  be  overwritten. 


Telephone:  (902)  492-4780 


FaedmUe:  (902)  492-4781 


Satlantic  Inc 


Richmond  Terminal,  Pier  9 
3295  Barrington  Street 
Halifax,  Hova  Scotia 
Canada  B3K  5X8 


5.2  Software 

A  disk  is  provided  with  each  system.  The  disk  contains  the  following  files: 

0CM_PR03.C  -  source  code  for  processing  program 
0CM_PR04.C  -  source  code  for  short  processing  program 
04252A.CAL  -  calibration  file  for  processing  data  from  OCM  4252 
[EXAMPLES]  -  examples  directory 

0CM_PR03.EXE  -  executable 

0CM_PR04.EXE  -  executable 

SEP05.ARG  -  sample  multiplatform  data  file 

INFILES. LST  -  sample  batch  processing  file 

04742A.CAL  -  sample  calibration  file 

04745A,CAL  -  sample  calibration  file 

04742  05. DAT  -  SEP05  processed  using  0CM_PR04 

04745~05.DAT  -  SEP05  processed  using  0CM_PR04 

The  file  SEP05.ARG  was  collected  during  our  JGOFS  aircraft  mission  in  1992, 
platform  04742  was  moored  in  Monterey  Bay,  California.  Platforms  04744  and  04745 
were  air-launched  over  the  equatorial  Pacific. 


Telephone:  (902)  492-^780  Facdmile:  (902)  492-4781 


Satlantic  Inc 


Richmond  Terminal,  Pier  9 
3295  Barrington  Street 
Halifax,  Hova  Scotia 
Canada  B3K  5X8 


6.0  General  Information 
6.1  Anti-Biofoulant: 

Toxic  anti'biofoulants  are  illegal  in  Canada,  thus  there  is  no  anti-biofoulant  on  the 
sensor.  However,  we  do  recommend  that  one  be  used,  if  it  is  legal  in  your  region.  The 
recommended  product  for  coating  optics  is  a  spray  can  of  Classic  Yacht  manufactured  by 
Philadelphia  Resins  which  is  sold  in  most  marine  stores  and  may  be  approved  by  the  EPA 
for  use  on  outboard  motors  in  some  areas.  Note  that  the  compound  contains  Tributytin 
Methacrylate  and  is  EXTREMELY  toxic.  The  warning  labels  should  be  read  carefully  and 
the  appropriate  respirator,  gloves,  goggles,  and  envirosuit  should  be  worn  Also  note  that 
some  batches  of  spray  cans  of  the  material  are  not  appropriate  for  use  on  optics.  You 
should  test  spray  on  a  clear  surface  and  let  dry.  if  it  appears  hazy  blue,  then  it  is  not 
appropriate.  Satlantic  is  not  responsible  for  the  results  when  using  such  compounds. 


6.2  Problems? 

If  you  have  any  problems  or  questions  about  the  instrument,  call  Satlantic  at 
(902)  492-4780  between  9am  and  5pm  Atlantic  Time,  or  FAX  us  at  (902)  492-4781. 


Telephone:  (902)  492-4780 


Facsimile:  (902)  492-4781 


Satlantic  Inc 


Richmond  Terminal,  Pier  9 
3295  Barrington  Street 
Halifax,  Nova  Scotia 
Canada  B3K  5X8 


APPENDIX  I 
Source  Code  Listings  for: 
0CM_PR03.C 
0CM_PR04.C 


Telephone:  (902)  492-4780 


Facsimile:  (902)  492-4781 


Program  0CM_PR03 . C 

This  program  ingests  datafiles  captured  from  the 
ARGOS  Processing  Center.  The  data  must  be  downloaded 
from  the  processing  center  in  TX  format  using  the  PRV 
command  (refer  to  ARGOS  Users  Manual  section  3.2.5. 

The  program  extracts  the  buoy  data  from  these  files 
and  produces  an  ASCII  file  with  calibrated  data  in 
physical  units. 


written  by:  Scott  McLean 
rlease  version  2.0 

ver  date:  Jan  15/93  for  SYMPHONY  formatted  files 

June  8/93  modified  for  ARGOS  TX  formatted  files 

distribution:  unlimited 

(c)  Satlantic  Inc.  1993 

Satlantic  Inc. 

3295  Barrington  St. 

Halifax,  NS 
CANADA  B3K  5X8 

telephone:  (902)  492-4780 
facsimile:  (902)  492-4781 


valid  for  buoy  data  from  05  Aug  92  -  date 


This  program  produces  output  files  which  are  very  wide 
(exceeding  256  characters),  if  this  is  too  large,  use  0CM_PR04.C 
which  does  not  include  min,  max  and  standard  deviation  values  (only 
the  averages)  and  produces  a  more  compact  output  file. 


This  program  reads  a  list  of  data  files  which  contain 
TX  formated  ARGOS  data  from  OCM  type  buoys.  Any  number  of  such 
files  can  be  read.  The  data  in  each  file  is  scanned  and  each 
data  transmission  will  be  processed,  all  extraneous  characters 
(for  example,  the  login)  are  ignored.  Only  one  platform  can  be 
processed  at  a  time,  if  you  have  three  platforms,  you  must  run 
program  three  times. 


Required  inputs: 

1)  INFILES. LST 

This  file  contains  a  list  of  all  files  that  are  to  be 
processed.  File  names  must  not  have  exetnsion  in  this  list. 

The  program  assumes  all  raw  data  files  have  the  .ARG  extension. 

2)  .ARG  files 

These  files  contain  the  data  obtained  on  various  days  from 


ms  =»»s  :|*s  =!►  =»ls 


the  ARGOS  processing  center  in  TX  format.  All  files  must  have 
the  .ARG  extension. 


3)  .CAL  files 

These  files  contain  the  calibration  data  for  each  platform. 

The  first  5  characters  of  this  file  name  are  the  platform  ARGOS 
ID  number.  The  program  will  prompt  you  for  the  calibration  file 
name.  It  will  remove  the  first  five  characters  from  the 
calibration  file  name  and  then  search  all  data  files  in  INFILES. LST 
for  data  hits  containing  this  ID. 

4)  output  file 

This  is  an  ASCII  output  file  containing  the  data  for  the  platfo 
requested.  The  program  will  prompt  you  for  a  file  name.  If  the  file 
already  exists,  it  will  be  overwritten. 


*/ 


include  <stdio.h> 
include  <conio.h> 
include  <math.h> 
include  <string.h> 
include  <stdlib.h> 

define  num  chans  7 
define  COl^ENT 
Sdefine  IRRAD  7 


/*  FUNCTION  PROTOTYPES  */ 

double  strtohex(char  hexchar) ; 

double  hextod(char  string[3]); 

double  cal_rad (double  counts,  int  channel) ; 

double  cal_irrad (double  counts) ; 

double  cal_temp( double  counts) ; 

double  explO (double  arg) ; 

double  stand_dev (double  sum_squares,  double  mean) ; 

double  get_tick (double  tick) ; 

double  get_yolts (double  tick); 

int  get_calibration_vals() ; 

int  hit(char  line[80],  char  platform[6] ) ; 


/*  GLOBAL  VARIABLES  */ 

double  scale(8] , offset [8] ;  /*  calibrations  */ 

FILE  *CAL_data; 
char  platfonB[6]; 


/*  FUNCTIONS  */ 


r 

I 


FUNCTION  STRTOHEX  converts  a  character  '0' 
lOuble  strtohex(char  hexchar) 
double  hex; 


.'F'  into  decimal  0..15  */ 


■  switch 

1 

(  hexchar  ) 

{ 

■ 

case 

/  i 

hex 

s 

0 

break; 

1 

case 

'0' 

hex 

= 

0 

break ; 

case 

'1' 

hex 

s: 

1 

break; 

case 

'2' 

hex 

= 

2 

break; 

1 

case 

'2' 

hex 

s= 

3 

break; 

■ 

case 

'4' 

hex 

= 

4 

break; 

case 

hex 

= 

5 

break; 

■ 

case 

'6' 

hex 

= 

6 

break ; 

1 

case 

/?/ 

hex 

7 

break; 

case 

,Q, 

hex 

= 

8 

break; 

case 

/g/ 

hex 

= 

9 

break; 

1 

case 

'k' 

hex 

2= 

10;  break; 

■ 

case 

hex 

11;  break; 

case 

'C' 

hex 

= 

12;  break; 

1 

case 

hex 

13;  break; 

case 

'E' 

hex 

=5 

14;  break; 

mm  \ 

case 

hex 

3S 

15;  break; 

1 

1  return 

hex; 

END 

FUNCTION 

STRTOHEX 

*/ 

/*  FUNCTION  HEXTOD  converts  a  string  '00'..'FF'  to  decimal  0..255  */ 

double  hextod(char  string[3]) 

{ 

double  hex; 

hex  «  strtohex ( string [ 0 ] )  *  16  +  strtohex(string[l] ) ; 
return  hex; 

/*  END  FUNCTION  HEXTOD  */ 

) 

/*  FUNCTION  GET_TICK  computes  time  since  last  transmit  buffer  update  */ 


double  get_tick (double  tick) 

{ 


/*  get  time  since  last  transmit  buffer  update  in  minutes  */ 

double  mins; 
unsigned  char  test; 

-est  =  tick; 

t.est  =  test  &  0x3F;  /*  mask  off  battery  voltage  in  upper  two  bits  */ 
mins  *=  test; 

mins  •=  mins  *  1;  /*  each  tick  is  1  minute  */ 
return  mins; 


*  END  FUNCTION  GET_TICK  */ 


FUNCTION  GET_VOLTS  calculates  current  battery  voltage  */ 


:  ouble  get_volts (double  tick) 


/*  battery  voltage  in  upper  two  bits,  each  bit  represents  1.5V  above 
10.5V  */ 


double  volts; 
msigned  char  test; 

test  «  tick; 
test  =  test  »  6; 
volts  =  test; 

volts  =  volts  *  1.5  +  10.5; 
return  volts; 


/*  END  FUNCTION  GET_VOLTS  */ 


/*  FUNCTION  CAL__RAD  calculates  upwelling  radiances  from  the  calibration 

values,  include  the  immersion  effects  */ 


double  cal_rad (double  counts,  int  channel) 

{ 

double  rad; 

double  immersiont]  =  {1.719,  1.720,  1.728,  1.732,  1.734,  1.741,  1.747); 


channel  *  channel  -  1;  /*  channel  o  =  Lu683  */ 

rad  »  scale [channel]  *  (counts  -  offset [channel ] )  *  immers ion [ channel ] ; 


^return  rad; 

/*  END  FUNCTION  CAL_RAD  */ 

I 

FUNCTION  CAL_IRRAD  calculates  downwelling  irradiance  from  mast  sensor 
*  calibration  values  */ 

liouble  cal_irrad (double  counts) 

■( 

I  ■  double  irrad; 

_  irrad  *  scale[IRRAD]  *  (counts  -  offset [ IRRAD] )  ;  /*  in  uW/cm''2/nin  */ 

■  return  irrad; 

1/*  END  FUNCTION  CAL_IRRAD  */ 

m) 

B/*  FUNCTION  CAL__TEMP  calculates  the  buoy  hull  temperature  */ 

■  double  cal^temp (double  counts) 

t 

.  double  temp ; 

■  temp  =  counts  *  0.16  -  5.0;  /*  degrees  C  */ 

return  temp; 

I  /*  END  FUNCTION  CAL_TEMP  */ 

■  } 

B  /*  FUNCTION  EXPIO  calculates  lO^X  */ 

B  double  explO (double  arg) 

{  double  val; 

B  val  «  exp  (  arg  *  log(lO)  ); 

_  return  val ; 


/*  END  FUNCTION  EXPIO  */ 


'*  FUNCTION  STAND_DEV  computes  the  encoded  standard  deviation  */ 


Ipiouble  stand_dev (double  suin_squares,  double  mean) 
/*  mean  should  be  in  volts  */ 


double  dev,  var; 
double  n  »  40.0; 
double  scale_fact  =  36.4; 

mean=“  3.825  /  255  *  mean; 
sum_squares  =  sum_squares  /  scale_fact; 

sum_squares  =  explO (sum_squares)  /  lOOOO;  /*  sum  of  squares  */ 
var  *  (  sum_squares  -  (mean  *  mean)  *  n  )  /  (n-l) ;  /*  variance  */ 

if  (var<0) 
dev=0 ; 
else 

dev  =  sqrt(var) ;  /*  standard  deviation  in  volts  */ 

dev  *  dev  *  255.0  /  3.825;  /*  convert  back  to  counts  for  cal_rad  */ 
return  dev; 


/*  END  FUNCTION  STAND_DEV  */ 


/*  PROCEDURE  PR0CESS_DATA_1  processes  line  1  of  a  data  hit  in  TX  format  */ 

void  process_data__l(char  line[803  , 

double  *iat, 
double  *lon, 
double  *class, 
double  *day, 
double  *pass_tiroe) 


char  lat__str  [  22  ]  ; 

char  lon__str  [8  ]=»•'”; 
char  class__str  [8 
char  day__str  [  8  ]  *”  ” ; 
char  pass_t  ime__str  [  8  ]  ; 

char  lat_hem; 
char  lon_hem; 
char  *end  str; 


/*  extract  parameters  from  input  string  */ 

sscanf  (line,  "%*8c%6c%c%*c%7c%c"  ,  lat__str , &lat_hem,  lon__str ,  &lon_hem)  ; 


Iscanf  (line, "%*27c%7s'' , class_str)  ; 
scanf  (line, ''%*40c%3s“ ,day_str)  ; 
sscanf  (line, "%*44c%4s” ,pass_tiine_str)  ;  /**/ 


convert  strings  to  numerics  */ 

|*lat  =  strtod(lat_str,  &end_str)  ; 

*lon  =  strtod(lon_str,  &end_str) ; 

♦class  =  strtod(class_str,  &end_str) ; 

|*day  =  strtod(day_str,  fiiend_str)  ; 

*pass_tiine  =  strtod(pass_time_str,  &end_str)  ; 

|*lat  *=  2*(lat_hem  «  'N')-l; 

♦Ion  ♦=  2^(lon~hem  *=  'E')--l; 

■*  END  procedure  PR0CESS_DATA_1  */ 


PROCEDURE  PR0CESS_DATA_2  processes  line  2  of  a  data  hit  in  TX  format  ♦/ 

void  process__data_2  (char  line [80], 

■double  mins(], 

double  maxs(], 
double  aves[], 
double  stddevs[])  /♦*/ 

{  char  chl_str[]  =  ”12345678"; 

I  char  ch2_str[]  =  "12345678"; 
char  ch3_str[]  =  "12345678"; 
char  ch4_str[]  =  "12345678"; 

■  char  inin_str[]  -  "  ",  max_str(]  =  "  ", 

ave_str[]  *  "  ",  stddev_str[]  =  "  "; 

char  ♦end  str; 


sscanf (line, "%*9c%8c%^5c%8c%*5c%8c%*5c%8c" , 

I  chl_str,  ch2_str,  ch3_str,  ch4_str) ; 

/♦  process  channel  1  Lu683  */ 

H  sscanf  (chl_str ,  "%2c%2c%2c%2c" ,  ave_str,roin__str , max__str ,  stddev_str) ; 

mins[l]  “  hextod(iain__str)  ; 

|maxs[l]  =  hextod (ttax_str ) ; 
aves(l]  *  hextod(ave_str)  ; 
stddevsCl]  *  hextod(stddev_str) ; 


/♦  process  channel  2  Lu670  ♦/ 


I  sscanf  ( ch2_str , "  %2c%2c%2c%2c" ,  ave_str ,  min_str ,  max__str ,  stddev_str )  ; 


rains[2]  »  hextod(min_str) ; 
iBaxs[2]  »  hextod(iaax_str)  ; 
aves{2j  »  hextod(ave_str) ; 
stddevs[2]  s=  hextod(stddev_str)  ; 

/*  process  channel  3  Lu555  */ 

sscanf (ch3_str,  *'%2c%2c%2c%2c'' , ave_str,inin_str,max_str, stddev_str)  ; 

mins [3]  =  hextod(min_str) ; 
maxs[3i  =  hextod(max_str) ; 
aves[3j  =  hextod(ave_str) ; 
stddevs[3]  ®  hextod(stddev_str) ; 

/*  process  channel  4  LuSlO  */ 

sscanf  (ch4_str,  "%2c%2c%2c%2c'* ,  ave_str,inin_str ,inax_str ,  stddev_str)  ; 

mins [4]  =  hextod(min_str) ; 
maxs[4i  =  hextod(roax_str) ; 
aves[4]  «  hextod(ave_str) ; 
stddevs[4]  *  hextod(stddev_str) ; 


/*  END  procedure  PR0CESS_DATA_2  */ 

) 


/*  PROCEDURE  PR0CESS_DATA_3  processes  line  3  of  a  data  hit  in  TX  format  */ 


void  proc€ss_data_3 


(char  line[80] , 

double  mins[], 
double  maxs [ ] , 
double  aves(3, 
double  stddevs [ ] , 
double  *irrad_ave, 
double  *irrad_std) 


/**/ 


{  char 
char 
char 
char 
char 

char 


chl_str[ ] 
ch2~str[ ] 
ch3~str[] 
ch4~str[ ] 
min~strC ] 
ave_str ( ] 
*end  str; 


»  "12345678”; 

=  "12345678"; 

=  "12345678"; 

*  "1234"; 

*  "  ",  max_str[3  ®  ” 

=  "  ",  stddev_str[3 


II 


II  • 


sscanf (line, "%*9c%8c%*5c%8c%*5c%8c%*9c%4c", 
chl^str,  ch2_str,  ch3_str,  ch4_str) ; 

/*  process  channel  5  Lu490  */ 

sscanf  (chl_str,  "%2c%2c%2c%2c" ,  ave_str ,min__str ,max__str ,  stddev__str)  ; 
min8[53  ■  hextod(min_str)  ; 


iBaxs[51  =  hextod(inax_str)  ; 
av€s{5]  =  hextod(ave_str) ; 
stddevs[5]  =  hextod(stddev_str) ; 

/*  process  channel  6  Lu443  */ 

sscanf  (ch2_str,  ''%2c%2c%2c%2c" ,  ave_str,inin_str ,  max_str,  stddev_str)  ; 

mins [6]  =  hextod(min_str) ; 
maxs[6]  =  hextod(max__str)  ; 
aves ( 6 ]  *  hextod ( ave_str ) ; 
stddevs[6]  «  hextod (stddev_str) ; 

/*  process  channel  7  Lu412  */ 

sscanf  {ch3_str,  ''%2c%2c%2c%2c" ,  ave_str,inin_str,  inax_str ,  stddev_str)  ; 

mins [7]  =  hextod (min^str) ; 
maxs[7]  =  hextod (max^str) ; 
aves [7]  s=  hextod (ave_str)  ; 
stddevs[7]  =  hextod (stddev_str) ; 

/*  process  irradiance  Ed490  */ 

sscanf  (ch4_str,  "%2c%2c" , ave__str , stddev^str)  ; 

*irrad_ave  =  hextod  (ave__str)  ; 

*irrad~std  =  hextod (stddev_str) ; 


/*  END  procedure  PR0CESS_DATA_3  */ 

) 

/*  PROCEDURE  PR0CESS_DATA_4  processes  line  4  cf  a  data  hit  in  TX  format  */ 


void  process__data_4  (char  line [80], 

double  *temp, 
double  *tick)  /**/ 


{  char  chl_str[]  =  ''12”; 
char  ch2_str[]  =  ”12''; 
char  *end  str; 


sscanf  (line, ''%*15c%2c%*llc%2c'',chl_str,  ch2_str)  ; 


/*  process  temperature  */ 

*temp  »  hextod (chl_str) ; 
♦tick  *  hextod (ch2_str) ; 


/*  END  procedure  PR0CESS_DATA_4  */ 


f  FUNCTION  LOCATION_CLASS  reports  accuracy  of  the  platform  location 

in  meters  from  the  ARGO 

l^auble  location_class (double  class) 

( 

I  double  accuracy; 

accuracy  =0.0; 

I  if  (class  ==  1.0) 

accuracy  =  1000.0; 
if  (class  ==  2.0) 

accuracy  =  350.0; 

I  if  (class  ==  3.0) 

accuracy  =  150.0; 
return  accuracy; 

END  FUNCTION  L0CATI0N_CLASS  */ 


FUNCTION  DEC_DATE  computes  the  decimal  day  from  the  day,  hour,  minute  */ 
louble  dec__date (double  day,  double  time) 


IP  double  hour,  minutes,  dec_day; 

hour  =  floor(time/100.0) ; 

■  minutes  =  time  -  hour  *  100; 

dec__day  =  day  +  hour/24.0  +  minutes/24.0/60.0; 
return  dec_day; 


■/*  END  FUNCTION  DEC_DATE  */ 


|/*  PROCEDURE  PROCESS_^HIT  uses  data  extracted  by  the  PROCESS_DATA_  procedures, 

applies  calibration  data  and  wr 

_  to  the  output  file  */ 


void  process_hit(FIIjE  *PROCESS__data, 

I”  double  day, 

double  time, 
double  lat, 

■double  Ion, 

doxible  class, 
double  aves [ ] , 


I  double  inins[], 

double  inaxs[], 
double  stddevs[], 

I  double  irrad_ave, 

double  irrad_stddev , 
double  temp, 

I  double  tick) 

7*  process  counts  and  produce  physical  units  for  each  parameter  */ 


double  accuracy; 

lint  channel, i; 

jdouble  std_dev ,  battery  ; 

double  ave,  min,  max,tmp,  minutes; 

.double  date; 

•accuracy*  location_class (class) ; 
date*dec_date(day , time) ; 


fprintf (PROCESS_data, "%s  %9.4lf  %6.31f  %6.3lf  %5. Olf" , platform, 

date , lat , Ion , accuracy) ; 

/*  apply  calibrations  to  radiance  channels  */ 

for  (channel*!;  channel  <  8;  channel++) 

I 

i=channel ; 
tmp*aves ( i ] ; 

ave  *  cal_rad(tmp,  channel) ; 
tmp=mins[i];  min  =  cal_rad(tmp,  channel); 
tmp=maxs[i];  max  =  cal_rad(tmp,  channel); 
tmp=stddevs [ i ] ; 

std__dev  =  stand__dev(tmp,  av€s(i]); 
std~dev  *  cal_rad(std_dev,  channel) ; 


fprintf (PROCESS_data, ”  %6.4lf  %6.4lf  %6.4lf  %6.4lf",  ave,  min,  max, 
~std__dev)  ; 

} 

/*  apply  calibration  to  irradiance  data  */ 

ave  *  cal_irrad(irrad_ave) ; 

std__dev  »  stand_dev(irrad_stddev,  irrad__ave) ; 

std__dev«  cal__irrad(std__dev) ; 

temp  »  cal_temp(temp) ; 
minutes  «  get_tick(tick) ; 
battery  «  get“volts(tick) ; 

fprintf (PROCESS_data,"  %5.11f  %5.11f  %4.11f  %4.11f  %4.11f\n'',ave,std_dev. 


temp, battery, minutes) ; 


i 


*  END  PROCEDURE  PROCESS_HIT 


*/ 


PROCEDURE  GET  CALIBRATION_VALS 


int  get_calibration_vals() 


int  error , i ; 
char  line [ 80 ]=••'•; 


reads  in  calibration  data  */ 


char  first  char; 


i=0; 

error=0 ; 

while  ( ! feof (CAL_data) ) 

{ 

line[0]  =  0; 

f gets  ( 1  ine  ,79,  CAL__data )  ; 

<}if  defined  DEBUG 
printf  ("%s*' , line)  ; 

#endi£ 


if  (strlen(line)  ==  o  |  1  lineCO]  ==  COMMENT)  continue; 
sscanf (line, "lie  %lf" ,scale+i,offset+i) ; 

#if  defined  DEBUG 

printf  (''%le  %lf\n", scale [i]  ,offset(i) )  ; 

#endif 

i++; 

if  ((i>8)  &&  (line(0]  "\n")) 

^printf (”  ERROR  in  calibration  file  -  too  many  cal  vals\n") ; 
error=-l; 

} 

) 

if  (i<8) 

^  printf ('•  ERROR  in  calibration  file  -  too  few  cal  vals\n»); 
error=-l; 

) 

fcIose(CAL_data) ; 
retuim  error; 


^  /*  END  PROCEDURE  GET_CALIBRATION_VALS  */  j 

j 


f*  PROCEDURE  HIT  */ 


int  hit(char  line[80], 

I  char  platfonn[6]) 


ll^har  checkt9]='"' ; 

if  (strstr (line, 

■  sscanf (line, "%1C 

I  if  ( strstr (checJ 
return  0; 


if  (strstr (line, platform)  ==  NULL)  return  0;  /*  if  ID  not  in  line  skip  */ 

sscanf  (line,  ”%10s'' ,  check)  ;  /*  check  first  10  chars  for  ID  only  */ 


if  (strstr (check,  platform)  ==  NULL) 
return  0; 
else 

return  1 ; 


I/**  PROCEDURE  MAIN  **/ 


I/** 

|main ( } 


I\ 

FILE  *ARGOS_data;  /*  current  ARGOS  file  */ 

FILE  *PROCESS_data ;  /*  output  file  */ 

FILE  *PROCESS_LIST;  /*  list  of  files  to  process  */ 

I  char  file_list[]  *  "infiles.  1st'' ; 
char  infile_ext[)  =  ".ARG"; 

■char  inf ile_name [ 13 ) ; 
char  outfile_naine[13]  ,* 
char  cal_f  ile__name  ( 13  ]  ; 

■char  line[80]  = 
int  line_no  *  0; 
int  data__hits  =  0; 

■double  lat, Ion, class, day, time; 

double  mins [8] ,  maxsCS],  aves[8],  stddevs[8]; 
double  irrad_ave,  irrad_std; 

■double  temp; 

double  tick_time; 
int  i; 


/*  read  in  list  of  files  to  .ARG  files  to  process  */ 
if  ( (PROCESS_LIST  »  fopen(file_list,"r"))  =  NULL) 

^printf("  File  INFILES. LST  must  exist  in  current  directory\n") ; 


printf("  and  must  contain  a  list  of  files  to  process i\n") ; 
return  -l; 

) 


/*  get  output  file  name  for  calibrated  data  */ 

puts{"  ENTER  output  file  name:"); 
scanf ("%13s",outfile_name) ; 

if  ( (PROCESS_data  =  fopen(outfile_name, "w") )  ==  NULL) 
( 

printfC  ERROR  creating  %s  \n'' ,outfile_na;ae)  ; 
return  -1; 

) 


/*  read  in  calibration  file  for  current  platform  */ 

putsC  ENTER  platform  calibration  file  name:"); 
scanf ("%13s" ,cal_file_name) ; 

if  ((CAL  data  =  fopen (cal_f ile_name^ "r") )  ==  NULL) 

{ 

printfC  ERROR  calibration  file  %s  not  found\n" ,  cal_file_name)  ; 
return  -1; 

) 


if  (get_calibration_vals 0  *=  -1) 

{ 

printf("  Program  Terminated  \n") ; 
return  -1;  /*  read  cal  file  */ 

) 


/*  extract  platform  name  from  first  5  characters  of  cal  file  name  */ 
sscanf  (cal_file_name,  "%5s‘' , platform) ; 


/*  put  text  header  on  output  file  */ 

fprintf (PROCESS_data , " 
fprintf (PROCESS_data, " 
fprintf ( PROCESS_data , " 


ID 

DAY 

LAT 

LON 

Lu683ave 

min 

max 

std" ) 

Lu670ave 

min 

max 

std" ) 

Lu555ave 

min 

max 

std") 

LuSlOave 

min 

max 

std") 

Lu490ave 

min 

max 

std") 

Lu443ave 

min 

max 

std" ) 

Lu412ave 

min 

max 

std") 

ACCUR") ; 


fprintf ( PROCESS_data , " 
fprintf ( PROCESS_data , ” 
fprintf ( PROCESS^data , " 

fprintf (PROCESS_data," _  _  _  _  , 

fprintf  (PROCESS_data,'*  ED490ave  std  Temp  Vbat  Tick\n 


’•) 


/*  scan  .ARG  files  and  extract  data  for  current  platform  */ 


while  (  lfeof(PROCESS__LlST)  ) 

{ 


line[0]®0; 


'  I  I  line(0]  -=  '\n^  )  continu 


r 

I 


infile_naine[0]=0; 
fgets(line,80,PROCESS_LIST) ; 

if  (  (line[0]  ==  ^\0')  \\  (line[0]  == 
sscanf  (line,  ”%s'' ,  inf ile_naine)  ; 


/*  scan  current  .ARG  file  for  current  platform  */ 


if  {(ARGOS_data  =  f open ( inf ile_name, "r”) )  !=  NULL) 

{ 

fgets(line,80,ARGOS_data) ;  /*  read  line  from  file  */ 

line  no++; 

while  (  Ifeof (ARGOS_data)  ) 

( 

/*  check  current  line  for  correct  plaform  ID  */ 

if  (hit(line, platform) )  /*  skip  line  if  not  data  hit  */ 

{ 

data_hits++; 

printf("  processing  hit  %3.0i  platform  %s\n", data_hits, platform) ; 
process__data_l (line,  &lat,  &lon,  &class,  &day,  &time)  ; 

f gets (line, 80, ARGOS_data) ; 

process_data_2 (line,  mins,  maxs,  aves,  stddevs) ; 
line_no++; 

f gets ( 1 ine ,80, ARGOS_data ) ; 

process_data_3 (line,  mins,  maxs,  aves,  stddevs,  &irrad_ave,  &irrad_st 
1 ine_no++ ; 

fgets(line,80,ARGOS_data) ; 
process_data_4 (line,  &terap,  &tick_time) ; 
line_no++; 

process_hit (PROCESS_data , day , time , lat , Ion , class , aves , mins , maxs , 
stddevs, irrad_ave,irrad_std, temp, tick  time) ; 

) 

fgets(line,80,ARGOS_data) ;  /*  read  next  line  */ 

line_no++; 

) 

line_no — ; 
fclose(ARGOS_data) ; 

} 

else 

printfC”  Error  file;  %s  not  found\n'' , inf ile_name)  ; 

}  /*  end  while  processing  */ 

fclose(PROCESS_LIST) ; 
fclose(PROCESS_data) ? 

printf("  %i  lines  processed\n",line_no) ; 
printf("  %i  data  hits  processed\n",data_hits) ; 
return  0; 


/***  end  KAIN  ***/ 


Program  0CM_PR04 . c 

This  program  ingests  datafiles  captured  from  the 
ARGOS  Processing  Center.  The  data  must  be  downloaded 
from  the  processing  center  in  TX  format  using  the  PRV 
command  (refer  to  ARGOS  Users  Manual  section  3.2.5. 

The  program  extracts  the  buoy  data  from  these  files 
and  produces  an  ASCII  file  with  calibrated  data  in 
physical  units. 


written  by:  Scott  McLean 
release  version  2.0 

ver  date:  Jan  15/93  for  SYMPHONY  formatted  files 

June  8/93  modified  for  ARGOS  TX  fonnatted  files 

distribution:  unlimited 

(c)  Satlantic  Inc.  1993 

Satlantic  Inc, 

3295  Barrington  St. 

Halifax,  NS 
CANADA  B3K  5X8 

telephone:  (902)  492-4780 
facsimile:  (902)  492-4781 


valid  for  buoy  data  from  05  Aug  92  -  date 

This  program  produces  output  files  which  only  contains  the 
average  values  to  produce  a  more  compact  output  file.  Use  0CM_PR03.C 
if  all  of  the  data  including  min,  max  and  stanadrd  deviation  is 
desired. 


This  program  reads  a  list  of  data  files  which  contain 
TX  formated  ARGOS  data  from  OCM  type  buoys.  Any  number  of  such 
files  can  be  read.  The  data  in  each  file  is  scanned  and  each 
data  transmission  will  be  processed,  all  extraneous  characters 
(for  example,  the  login)  are  ignored.  Only  one  platform  can  be 
processed  at  a  time,  if  you  have  three  platforms,  you  must  run 
program  three  times. 


Required  inputs: 

1)  INFILES. LST 

This  file  contains  a  list  of  all  files  that  are  to  be 
processed.  File  names  must  not  have  exetnsion  in  this  list. 

The  program  assumes  all  raw  data  files  have  the  .ARG  extension. 

2)  .ARG  files 

These  files  contain  the  data  obtained  on  various  days  from 


FUNCTION  STRTOHEX  converts  a  character  '0‘ 


F'  into  decimal  0..15  */ 


pble  strtohex (char  hexchar) 

[ouble  hex; 

witch  (  hexchar  ) 

( 


case 

/  / 

hex 

0;  break; 

case 

'0' 

hex 

= 

0 ;  break ; 

case 

'l' 

hex 

= 

1;  break; 

case 

'2' 

hex 

= 

2;  break; 

case 

'3' 

hex 

= 

3;  break; 

case 

'4' 

hex 

= 

4 ;  break ; 

case 

'5' 

hex 

5;  break; 

case 

'6' 

hex 

= 

6;  break; 

case 

'7' 

hex 

= 

7;  break; 

case 

'8  ' 

hex 

= 

8;  break; 

case 

'9' 

hex 

9;  break; 

case 

'A' 

hex 

= 

10;  break; 

case 

'B' 

hex 

= 

11;  break; 

case 

'C' 

hex 

12;  break; 

case 

'D' 

hex 

= 

13;  break; 

case 

'E' 

hex 

= 

14;  break; 

case 

\ 

'F' 

hex 

s 

15;  break; 

) 

return  hex; 

(*  END  FUNCTION  STRTOHEX  */ 


FUNCTION  HEXTOD  converts  a  string  '00'.. 'FF'  to  decimal  0..255  */ 


iouble  hextod(char  string[3]) 


double  hex; 

hex  =  strtohex (string[ 0] )  *  16  +  strtohex ( string [ 1 ]) ; 
return  hex; 


/*  END  FUNCTION  HEXTOD  */ 


FUNCTION  GET_TICK  computes  time  since  last  transmit  buffer  update  */ 


Iouble  get__tick  (double  tick) 


get  time  since  last  transmit  buffer  update  in  mimutes  */ 

double  mins; 
msigned  char  test; 

test  *  tick; 

test  =  test  &  0x3F;  /*  mask  off  battery  voltage  in  upper  two  bits  */ 
mins  *  test; 

mins  *  mins  *  1;  /*  each  tick  is  l  minute  */ 
return  mins; 


*  END  FUNCTION  GET_TICK  */ 


*  FUNCTION  GET_VOLTS  calculates  current  battery  voltage  */ 


’.ouble  get_volts  (double  tick) 


/*  battery  voltage  in  upper  two  bits,  each  bit  represents  1.5V  above 
10.5V  */ 


double  volts; 
unsigned  char  test; 

test  =  tick; 
test  =  test  >>  6; 
volts  *  test; 

volts  «  volts  *  1.5  +  10.5; 
return  volts; 


/*  END  FUNCTION  GET_VOLTS  */ 

) 

/*  FUNCTION  CAL_RAD  calculates  upwelling  radiances  from  the  calibration 

values,  include  the  immersion  effects  */ 


double  cal_rad (double  counts,  int  channel) 

I 

double  rad; 

double  immersion[]  =  (1.719,  1.720,  1.728,  1.732,  1.734,  1.741,  1.747) 
channel  «*  channel  -  1;  /*  channel  0  »=  Lu683  */ 


rad  -  scale [channel]  *  (counts  -  offset(channel] )  *  immersion [channel] 


return  rad; 


END  FUNCTION  CAL_RAD  */ 


FUNCTION  CAL_IRRAD  calculates  downwelling  irradiance  from  mast  sensor 
I  calibration  values  */ 


double  cal_irrad (double  counts) 


double  irrad; 


Hirrad  =  scale[IRRAD]  *  (counts  -  offset [ IRRAD] ) ;  /*  in  uW/cm''2/nm  */ 


return  irrad; 


/*  END  FUNCTION  CAL_IRRAD  */ 


"/*  FUNCTION  CAL_TEMP  calculates  the  buoy  hull  temperature  */ 


louble  cal_temp (double  counts) 


II 

double  temp; 

temp  *=  counts  *  0,16  -  5.0;  /*  degrees  C  */ 

I  return  temp; 

/*  END  FUNCTION  CAL_TEMP  */ 

I, 

FUNCTION  EXPIO  calculates  10 */ 

I  double  explO (double  arg) 

I  {  double  val; 

val  *  exp  (  arg  *  log (10)  ) ; 

I  return  val ; 

|/*  END  FUNCTION  EXPIO  */ 


/*  FUNCTION  EXPIO  calculates  10 ^X  */ 


double  explO (double  arg) 


FUNCTION  STAND_DEV  computes  the  encoded  standard  deviation  */ 


louble  stand_dev  (double  sui!i_squares,  double  mean) 


mean  should  be  in  volts  */ 


double  dev,  var; 

■double  n  =  40.0; 
double  scale_fact  =  36.4; 

■mean®  3.825  /  255  *  mean; 

sum_squares  =  suni_squares  /  scale_fact; 

sum^squares  ®  explO (suiii_squares)  /  lOOOO;  /*  sum  of  squares  */ 
var  =  (  sum_squares  -  (mean  *  mean)  *  n  )  /  (n-1) ;  /*  variance  */ 

I  if  (var<0) 
dev=0 ; 
else 

H  dev  =  sqrt(var) ;  /*  standard  deviation  in  volts  */ 

dev  =  dev  *  255.0  /  3.825;  /*  convert  back  to  counts  for  cal_rad  */ 
_  return  dev; 


'*  END  FUNCTION  STAND_DEV  */ 


PROCEDURE  PR0CESS_DATA_1  processes  line  1  of  a  data  hit  in  TX  format  */ 

void  process_data__l(char  line[80], 

I  double  *lat, 

I  double  *lon, 

double  *class, 
double  *day, 
double  *pass_tiroe) 

{ 

char  lat_str[  22  ]=•"•; 
char  lon_str  [  8  ] ; 
char  class_str[8  ]='"•; 

I  char  day__str[8]='"* ; 

I  char  pass_time__str[8]='"' ; 

char  lat_hem ; 

,  char  lon^hem; 
char  *end  str; 


I  /*  extract  parameters  from  input  string  */ 

sscanf  (line, ''%*8c%6c%c%*c%7c%c" , lat_str, Silat^hem, lon_str , &lon_hem)  ; 


sscanf  (line, ''%*27c%7s'* , class_str)  ; 

sscanf (line, "%*40c%3s" ,day_str) ; 

sscanf  (line, ''%*44c%4s" ,  pass_tiine_str)  ;  /**/ 


/*  convert  strings  to  numerics  */ 

*lat  =  strtod(lat_str,  &end_str) ; 

*lon  =  strtod(lon_str ,  &end_str) ; 

♦class  »»  strtod(class_str,  &end_str) ; 

♦day  *  strtod(day_str,  &end_str)  ; 

♦pass^time  *  strtod(pass_tiine_str,  &end_str)  ; 

♦lat  ♦=  2^(lat_hem  «  'N')-l; 

♦Ion  ♦*  2^(lon_hem  ==  'E')-l; 


[/♦  END  procedure  PR0CESS_DATA_1  */ 


/♦  PROCEDURE  PR0CESS_DATA_2  processes  line  2  of  a  data  hit  in  TX  format  */ 


void  process_data__2  (char  line [80], 

double  mins [ ] , 
double  maxs[], 
double  aves[], 
double  stddevs[])  /**/ 


{  char 
char 
char 
char 
char 

char 


chl_str[ ] 
ch2_str[ ] 
ch3_str[ ] 
ch4__str[  ] 
min_str[ ] 
ave_str [ ] 
♦end  str; 


=  ”12345678''; 

=  ”12345678"; 

=  ”12345678"; 

«  "12345678"; 

=  ”  ",  max_str[]  =  " 

=  "  ",  stddev_str [ ] 


It 


It 


sscanf (line, "%*9c%8c%*5c%8c%^5c%8c%^5c%8c" , 
chl_str,  ch2_str,  ch3_str,  ch4_str) ; 

/♦  process  channel  i  Lu683  ♦/ 

sscanf (chl_str , "%2c%2c%2c%2c" , ave_str ,inin_str , max^str , stddev_str) ; 

minsfl]  “  hextod(iDin_str) ; 
maxs[l]  *  hextod(max_str) ; 
aves[l]  =  hextod(ave_str); 
stddevs[l]  »  hextod(stddev_str) ; 


/*  process  channel  2  Lu670  ♦/ 

sscanf  (ch2__str, **%2c%2c%2c%2c",ave_str,min_str,inax_str, stddev_str) ; 


|inins(2]  «  hextod{min_str)  ; 
maxs[2]  «  hexto<i(inax_str)  ; 
aves[2]  »  hextod(ave_str) ; 

||Stddevs[2]  «  hextod(stddev_str)  ; 

/*  process  channel  3  Lu555  */ 

m  sscanf  (ch3_str^  '•%2c%2c%2c%2c" ,  ave_str, inin_str ,inax_str ,  stddev_str)  ; 

mins[3]  =  hextod(inin_str)  ; 

|maxs[3]  “  hextod(inax_str)  ; 
aves [ 3 ]  »  hextod ( ave_str ) ; 
stddevs[3]  =  hextod (stddev_str) ; 

I  /*  process  channel  4  Lu510  */ 

sscanf  (ch4_str ,  ”%2c%2c%2c%2c" , ave_str,min_str ,inax_str , stddev_str)  ; 

I  mins [4]  *=  hextod  (min_str)  ; 
niaxs[4]  *  hextod (niax_str)  ; 

I  aves [4]  =  hextod (ave_str) ; 

stddevs[4]  »  hextod (stddev_str) ; 


I'*  END  procedure  PR0CESS_DATA_2  */ 


/*  PROCEDURE  PR0CESS_DATA_3  processes  line  3  of  a  data  hit  in  TX  format  */ 


void  process_data_ 


3 (char  line [80], 

double  mins[3, 
double  maxs[], 
double  aves [ ] , 
double  stddevs [ ] , 
double  *irrad_ave, 
double  *irrad_std)  /**/ 


{  char  chl_str[]  » 
char  ch2_str(]  - 
char  ch3_strn  * 
char  ch4_strt]  ’ 
char  min_str(]  = 
ave__str  [  ] 
char  *end  str; 


•'12345678"; 

"12345678"; 

"12345678"; 

"1234"; 

"  ",  max_str[]  =  «  ", 

=  "  ",  stddev_str[]  ®  "  ", 


sscanf (line, "%*9c%8c%*5c%8c%*5c%8c%*9c%4c" , 
chl__str,  ch2_str,  ch3_str,  ch4__str) ; 

/*  process  channel  5  Lu490  V 

sscanf  (chl_str,  "%2c%2c%2c%2c"  ,ave__str,min__str,max_str,  stddev_str)  ; 
nlns[5]  ■  hextod (inin_str) ; 


Btaxs[5]  »  hextod(inax_str)  ; 
aves[S]  *  hextod(ave_str) ; 
stddevs[5]  =  hextod(stddev_str) ; 

/*  process  channel  6  Lu443  */ 

sscanf  (ch2_str,  "%2c%2c%2c%2c" ,  ave_str,inin_str,max_str,  stddev_str)  ; 

mins [6]  =  hextod(min_str) ; 
maxs[6]  »  hextod(max_str) ; 
aves[61  =>  hextod(ave_str)  ; 
stddevs[6]  =  hextod(stddev_str) ; 

/*  process  channel  7  Lu412  */ 

sscanf  (ch3__str,  ”%2c%2c%2c%2c" ,  ave_str,min_str,roax_str,  stddev_str)  ; 

mins [7]  =s  hextod(min_str)  ; 
maxs[7]  =  hextod(max_str) ; 
aves[7]  =  hextod(ave_str) ; 
stddevs[7]  =  hextod(stddev_str) ; 

/*  process  irradiance  Ed490  */ 

sscanf  (ch4_str,  '’%2c%2c"  ,  ave_str,  stddev_str)  ; 

*irrad_ave  =  hextod(ave__str)  ; 

*irrad”std  =»  hextod(stddev_str)  ; 


/*  END  procedure  PR0CESS_DATA_3  */ 

) 

/*  PROCEDURE  PR0CESS_DATA_4  processes  line  4  of  a  data  hit  in  TX  format  */ 


void  process_data_4 (char  line[80] , 

double  *temp, 
double  *tick)  /**/ 


{  char  chl_str[]  ®  ”12"; 
char  ch2_str[]  =  "12  ■ 
char  *end  str; 


sscanf (line, "%*15c%2c%*llc%2c" ,chl_str,  ch2_str) ; 


/*  process  temperature  */ 

*temp  *  hextod(chl_str) ; 
♦tick  *»  hextod(ch2~str) ; 


/*  END  procedure  PR0CESS_DATA_4  */ 


I 

I 

FUNCTION  LOCATION_CLASS  reports  accuracy  of  the  platform  location 
I  in  meters  from  the  ARGO 

Ijdouble  location_class (double  class) 

{ 

H  double  accuracy ; 

*  accuracy  *0.0; 

I  if  (class  ==  1.0) 

accuracy  =  1000.0; 
if  (class  2.0) 

accuracy  =  350.0; 

I  if  (class  ==  3.0) 

accuracy  =  150.0; 
return  accuracy; 


/*  END  FUNCTION  LOCATION_CLASS  */ 


/*  FUNCTION  DEC_DATE  computes  the  decimal  day  from  the  day,  hour,  minute  */ 


double  dec__date (double  day,  double  time) 

{ 

double  hour,  minutes,  dec^day; 

hour  *  floor (time/100. 0) ; 
minutes  =  time  -  hour  *  100; 

dec_day  =  day  +  hour/24.0  +  roinutes/24 . 0/60. 0; 
return  dec_day; 


/*  END  FUNCTION  DEC_DATE  */ 

) 


/*  PROCEDURE  PROCESS_HIT  uses  data  extracted  by  the  PROCESS_DATA_  procedures, 

applies  calibration  data  and  wi 
to  the  output  file  */ 


void  process_hit(FILE  *PROCESS_data , 

double  day, 
double  time, 
double  lat, 
double  Ion, 
dovible  class, 
double  aves [ ] , 


double  mins[], 
double  roaxs[], 
double  stddevs[ ] , 
double  irrad_ave, 
double  irrad_stddev, 
double  temp, 
double  tick) 

/*  process  counts  and  produce  physical  units  for  each  parameter  */ 


double  accuracy; 

lint  channel, i; 

double  std__dev ,  battery ; 

double  ave,  min,  max,tmp,  minutes; 

double  date; 

accuracy*  location_class (class) ; 
date*dec_date (day, time) ; 


f print f (PROCESS^data, "%s  %9.41f  %6.31f  %6.31f  %5. Olf" , platform, 

date, lat, Ion, accuracy) ; 

/*  apply  calibrations  to  radiance  channels  */ 

for  (channel*!;  channel  <  8;  channel++) 

i 

i=channel ; 
tmp=aves [ i ] ; 

ave  =  cal_rad(tmp,  channel) ; 
tmp*roins[i];  min  =  cal_rad(tmp,  channel); 
tmp=maxs[i];  max  *  cal_rad(tmp,  channel); 
tmp*stddevs [ i ] ; 

std_dev  =  stand_dev(tmp,  aves[i]); 
std_dev  =  cal_rad(std_dev,  channel) ; 


fprintf  (PROCESS_data,"  %6.41f'',  ave)  ; 

) 

/*  apply  calibration  to  irradiance  data  */ 

ave  «  cal_irrad(irrad__ave) ; 

std__dev  *  stand__dev(irrad_stddev,irrad__ave) ; 

std_dev*  cal_irrad(std__dev) ; 

temp  «*  cal_temp(temp) ; 
minutes  «  get_tick{tick) ; 
battery  «  get“volts(tick) ; 

fprintf (PROCESS_data,"  %5.11f  %4.11f  %4.llf  %4.11f\n",ave, 

temp, battery# minutes) ; 


I 

I 

END  PROCEDURE  PROCESS_HIT  */ 


PROCEDURE  GET_CALIBRATION_VALS  reads  in  calibration  data  */ 


nt  get_calibration_vals() 


int  error,!; 
char  line[80]=""; 

char  first  char; 


1=0; 

error=0 ; 

while  ( ! feof (CAL_data) ) 

{ 

line[0]  =  0; 

f gets  ( 1  ine ,  7  9 ,  CAL__data )  ; 

#if  defined  DEBUG 
printf  (''%s"  ,  line)  ; 

#endif 


if  (strlen(line)  ==  o  | |  lineCO]  =-  COMMENT)  continue; 
sscanf (line, ”%le  %lf",scale+i,offset+i) ; 

#if  defined  DEBUG 

printf (”%le  %lf\n”,scale[i] ,offset[i]) ; 

#endif 

i++; 

if  ((i>8)  &&  {line[0]  1=  '\nM) 

printf (”  ERROR  in  calibration  file  -  too  many  cal  vals\n”) ; 
error=-l ; 

) 

) 

if  (i<8) 

^  printf ("  ERROR  in  calibration  file  -  too  few  cal  vals\n") ; 
error*-l ; 

) 

fclose(CAL_data) ; 
return  error; 


/*  END  ^PROCEDURE  GET_CALIBRATION_VALS  */ 


m*  PROCEDURE  HIT  */ 

Int  hit(char  line[80], 

char  platform [ 6 ] ) 

l( 

char  ch€ck[9]=""; 

I 

if  (strstL (line, platform)  =*  NULL)  return  0;  /*  if  platform  ID  not  in  line 

■  sscanf  (line, ''%10s'* ,  check)  ;  /*  check  first  10  chars  for  ID  only  */ 

if  (strstr (check,  platform)  ==  NULL) 

I  return  0; 

■  else 


return  1 ; 


END  PROCEDURE  PROCESS  HIT  */ 
|/**  PROCEDURE  MAIN  **/ 

HSiain  ( ) 


■FILE  *ARGOS_data;  /*  current  ARGOS  file  */ 

FILE  *PROCESS_<iata ;  /*  output  file  */ 

FIUE  *PROCESS_LIST;  /*  list  of  files  to  process  */ 

■char  file_list(]  =  " infiles. 1st" ; 
char  infile_ext[3  *  ".ARG”; 
char  infile~name[13] ; 

■  char  out  f i 1 e_naroe [13]; 
char  cal_file_name[13] ; 
char  line [80]  = 
int  line_no  «  0; 

I  int  data”hits  =  0; 

double  lat, Ion, class, day, time; 

■double  mins [8] ,  maxs[8],  aves[8],  stddevs[8]; 
double  irrad_ave,  irrad__std; 
double  temp; 

_  do\ible  tick  time; 

■  int  i; 


/*  read  in  list  of  files  to  .ARG  files  to  process  */ 
if  ( (PROCESS_LIST  «  fopen(file__list,«r'*))  NULL) 


I 

I 


{ 

printfC  File  INFILES. LST  roust  exist  in  current  directory\n" 
printf("  and  must  contain  a  list  of  files  to  process! \n")  ; 
return  -1; 

) 


|/*  get  output  file  name  for  calibrated  data  */ 


puts("  ENTER  output  file  name:")*' 
scant (”%13s",outfile_naroe) ; 

if  ( (PROCESS_data  =  fopen(outfile_name, '‘w**) )  ==  NULL) 
{ 

printfC  ERROR  creating  %s  \n'' ,outfile_name)  ; 
return  -1; 

) 


/*  read  in  calibration  file  for  current  platform  */ 

puts(’'  ENTER  platform  calibration  file  name:"); 
scant C%13s",cal_file_name)  ; 

if  ((CAL_data  *  fopen(cal  f ile_name, “r") )  ==  NULL) 

{ 

printfC  ERROR  calibration  file  %s  not  found\n" ,  cal_f  ile_name)  ; 
printf ( "  Program  Terminated  \n" ) ; 
return  -1; 

} 


if  (get_calibration  vals()  ==  -1) 

{ 

printf ("  Program  Terminated  \n")  ; 
return  -1;  /*  read  cal  file  */ 

} 


/*  extract  platform  name  from  first  5  characters  of  cal  file  name  */ 
sscanf (cal_file_name, "%5s", platform) ; 


/*  put  text  header  on  output  file  */ 


f printf  (PROCESS  data,"  ID  DAY  LAT  LON 

fprintf (PROCESS-data , "  Lu683 " ) ; 


f print f ( PROCESS  data , " 
fprintf (PROCESS“data, " 
fprintf (PROCESS_data , " 
f print f ( PROCESS  data , " 
fprintf (PROCASS”data, " 
fprintf (PROCESS”data, " 
fprintf (PROCESS_data , " 


LU670") ; 

Lu555")  ; 

LuSlO")  ; 

Lu490")  ; 

LU443")  ; 

Lu412")  ; 

ED490  Temp  Vbat  Tick\n") ; 


ACCUR" ) ; 


/*  scan  .ARG  files  and  extract  data  for  current  platform  */ 


Jhile  (  Ifeof  (PROCESS_LIST)  ) 
line[0]=0 ; 

I  inf iie_name [0]=0; 

fgets(line,80,PROCESS_LIST) ; 

if  (  (line[0]  ==  '\0')  11  (linetO]  ==  '  ')  11  line[0]  ==  '\n'  )  continu 
■  sscanf  (line,  "%s*’ ,  inf  ile__name)  ; 


/*  scan  current  .ARG  file  for  current  platform  */ 


if  ((ARGOS_data  »  fopen(infile_name, "r") )  !=  NULL) 

{ 

fgets(line, 80, ARGOS_data) ;  /*  read  line  from  file  */ 

line_no++; 

while  (  Ifeof (ARGOS  data)  ) 

{ 

/*  check  current  line  for  correct  plaform  ID  */ 

if  (hit (line, platform) )  /*  process  data  */ 

{ 

data_hits++; 

printfC’  processing  hit  %3.0i  platform  %s\n'' ,data_hits, platform)  ; 
process_data_l(line,  &lat,  &lon,  &class,  &day,  £<time); 

fgets(line,80,ARGOS_data) ; 

process_data_2 (line,  mins,  maxs,  aves,  stddevs) ; 

1 ine_no++ ; 

fgets(line,80,ARGOS_data) ; 

process_data_3 (line,  mins,  maxs,  aves,  stddevs,  &irrad_ave,  &irrad_st 
line_no++; 

f gets ( 1 ine ,80, ARGOS_data ) ; 
process_data_4 (lineT  Stemp,  &tick_time) ; 

line_no++; 

process__hit  (PROCESS_data ,  day ,  time ,  lat ,  Ion ,  class ,  aves ,  mins ,  maxs , 
stddevs, irrad_ave, irrad_std, temp, tick_tiroe) ; 

) 

fgets(line,80,ARGOS_data) ;  /*  read  next  line  */ 

line_no++; 

) 

line_no — ; 
f close (ARGOS_data) ; 

) 

else 

printf("  Error  file;  %s  not  found\n” , inf ile_name) ; 

)  /*  end  while  processing  */ 

fclose(PROCESS_LIST) ; 
fclose(PROCESS_data) ; 

printf("  %i  lines  processed\n",line_no) ; 
printf^^"  %i  data  hits  processed\n",data_hits) ; 
return  O; 
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APPENDLX  2 

Notes  on  TX  formatted  data  from  the  ARGOS  Users  Manual 


Telephone:  (902)  492-4780  Facsimile:  (902)  492-4781 


3.2.5  -  DOWNLOADING  RESULTS  IN  TX  FORMAT 


Command  PRV  lets  you  obtain  the  results,  in  the  TX 
format,  generated  by  one  or  more  platforms  in  a  program  or 
programs  of  which  you  are  the  owner.  The  results  for  the 
current  day  and  the  four  previous  days  can  be  provided. 

Use  command  PRV  to  select  the  results  you  wish  to 
download  by  entering  the  platform  numbers  and  message 
I  reception  dates. 


I 

i 

t 

I 

Specifies  the  type  of 
format  requirec 

Selected  program  number 

Provides  result  in  TX  format 

285/12  selects  results 
received  since  that  date, 
where  285  is  the  calendar 
date  and  12;00  the  time. 
Default  value;  current  day 
at  00:00. 

285/20  selects  the  results 
received  until  that  date, 
where  285  is  the  calendar 
date  and  20:00  the  time. 
Default  value:  current  day 
and  time. 

10000  is  the  chosen 


COMMANDPRV 

Enter  the  command  as  follows: 


platform  number. 


The  above  command  was  executed  on  September  13, 
1987  (calendar  day;  286).  The  selected  dates,  entered  in 
the  command  (285/12  -  285/20),  are  obviously  examples 
only  and  will  depend  on  the  date  of  your  interrogation. 
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USER  MANUAL 


Press  Carriage  Return;  the  following  results  appear; 


r 

ARGOS  READY 

/  PRV,10,TX,2a5/12~285/20.10000 
PROG  00010 


These  three  sets  of  results 
represent  the  most 
significant  messages 
received  on  the  285:h  day 
between  12;00  and  20  00 
The  structure  of  the  three 
sets  is  identical 


10000 

(2) 

43.543N  1.40 IE  3 

0.t0466E+4 

00 

28V1340Z-285/1336 
0.20000E+2  0  32000E+2 

10000 

43.544N  1.399E  3 

285/1521Z-285/1515 

(12) 

0.99875E*3 

32 

0.20000E+2  0.31000E*2 

\l 

lioooo 

I  (6) 

43.544N  1 .399E  3 

0.99927E+3 

32 

285/1 701Z-285/1S15 
0.20000E+2  0.30000E.^2 

ARGOS  READY 
/ 

V _ J 

Interpret  the  above  results  as  follows: 


Location  class 

Platform  longitude  in  degrees  and 
thousandths  of  a  degree 

Selected  program  number 

Selected  platform  number 


Calendar  day  and  time  (HHMM)  of 
last  location  in  UTC. 

Calendar  day  and  time  (HHMM)  of 
data  collection  in  UTC. 

Compression  index. 
Twelve  identical  messages 
were  received. 

Platform  latitude  in  degrees  and 
thousandths  of  a  degree 

Physical  value  for  sensor  l 
(998.75  hPa) 

Hexadecimal  value  for  sensor  2 
Physical  value  for  sensor  3  (20'C) 
Physical  value  for  sensor  4  (31  "C) 
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I 

I 


1  To  request  on-line  information  on  command  PRV, 
enter  "?,PRV'’  and  a  Carriage  Return. 

To  abort  the  current  command,  enter  "A". 

To  interrupt  the  current  command,  enter  "I"  (upper¬ 
case).  To  resume,  press  Carriage  Return. 

Command  PRV  can  be  shortened  to  "P". 

The  examples  below  show  the  options  available  under 
command  PRV. 

IPRV,10,TX„ 

Downloads  all  messages  processed  for  program  10 
received  on  the  current  day  since  00.00  UTC. 

IPRV,10,TX,  285, 

Downloads  all  messages  processed  for  program  10 
received  between  00.00  UTC  on  day  285  and  the  current 
day.  In  practice,  only  the  results  of  the  four  days  prior  to 
the  current  day  are  accessible. 

IP, 10, TX, 285118, 

Downloads  all  messages  processed  for  program  10 
received  between  day  285  at  18.00  and  the  current  day.  In 
practice,  only  the  results  of  the  four  days  prior  to  the 
current  day  are  accessible. 


IP„TX, 285-286, 

Downloads  all  messages  processed  (for  all  platforms  user 
has  access  to)  received  between  00.00  on  day  285  and 
24.00  on  day  286. 
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APPENDLX  3 

Using  0CM_PR04.C  an  example 

Includes:  sample  data  file  SEP05.ARG 

sample  calibration  file  04742A.CAL 
sample  calibration  file  04745 A.  CAL 
sample  batch  processing  file  INFILES. LST 
sample  of  0CM_PR04  executing  on  SEP05.ARG 
sample  output  4742_05.DAT 
sample  output  4745_05.DAT 


Telephone:  (902)  492-4780 


Facsimile:  (902)  492-4781 


\USR\SCOTT\EXPERIME .NTS\ JGOFS\BUOY\SEP05 . ARG 
Wednesday  July  28,  1993  07:28:04  pm 


Page:  1 


September  06,  1992  15:48  hr 
LOGIN  AT  249/1841  LAST  ACCESS  AT  248/1345  UTC 

ARGOS  READY 

/PRV, 895, TX, 248/ 13-249/ 18. 

SYNTAX  ERROR 

ARGOS  READY 

/PRV, 895, TX, 248/ 13-249/18, 

Prog  00895 

04742  36,748N  121.865W  0  248/1353Z-248/1352 

(  1)  00  35  21  10 

11  16  06  14 


04742 

36,748N  121.865W 

0 

(  1) 

00 

35 

11 

16 

7B 

CF 

04742 

36.746N  121.864W 

2 

(  1) 

13B 

145 

200066E 

7A 


200056A 

El 


04742 
<  1) 


36.745N  121.863W  3 

4030679  4020573 

191120AD  140E1AA7 

7A  CD 


04742  36.745N  121.863W  3 

(  1)  6050883  504067D 

FFFFFFFF  FFFFFFFF 

FI  C7 

04742  36.747N  121.867W  2 

(  1)  8060A8B  7050886 

2F2237C1  281F2EBC 

7B  DE 

04742  36,747N  121.866W  3 

(  1)  7040B88  6040980 

24163BB9  1E1327B3 

7B  C9 

04742  36.747N  121.866W  3 

(  1)  6030881  4020678 

1A0D22AE  140A1CA7 

7B  D4 

04742  36.746N  121.867W  2 

(  1)  100024D  100024E 

3000773  200056C 

7B  Cl 


04742  36.749N  121.864W 


248/1510Z-248/1506 
1000356  100035E 

2000464  268 


248/1649Z-248/1646 
C080F96  F0A139E 

120C17A3  12A3 


248/1824Z-248/1646 
110F14A2  167E003F 

FF76FFFF  FFFF 


248/2202Z-248/2158 
18121DAC  1D1622B2 

251E2EBA  21B6 


248/2340Z-248/2339 
130B21A5  170E27AB 

1B1225B0  lAAE 


249/0047Z-248/2339 
E06129B  100816A0 

130919A5  13A6 


249/0230Z-249/0226 
2000461  2000566 

2000569  370 


249/0410Z-249/0407 


r 
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(  1)  00 
FF55FFE7 
21 


34 

FFFFDEC3 

39 


22 

BFEF7BF9 


15 

EFC6 


04742 

36.747N  121.863W 

2 

(  1) 

00 

34 

OD 

14 

7B 

CO 

04742 

36.745N  121.867W 

2 

(  1) 

00 

34 

OC 

12 

7B 

ED 

04742 

36.750N  121.868W 

1 

(  1) 

00 

34 

OF 

14 

7B 

D5 

04742 

36.745N  121.862W 

2 

(  1) 

ID 

3A 

24B 

249 

7B 

DC 

04742 

36.755N  121.848W 

0 

(  1) 

3020572  302036B 

E09129C  E030F96 

C6 

38 

04744 

0.165S  140.958W 

2 

{  1) 

OA 

00 

27 

31 

BA 

D9 

04744 

0.165S  140.958W 

2 

(  1) 

OA 

00 

27 

38 

BA 

E7 

04744 

0.165S  140.958W 

2 

(  1) 

23E 

242 

A021795  A031796 

BA 

DC 

04744 

0.165S  140.958W 

2 

(  1) 

2010361  2010461 

2C1A38BF  30153DC2 

249/1023Z-249/1019 
21  OE 

00  OF 


249/1207Z-249/1202 
21  OF 

00  OE 


249/1342Z-249/134C 
IF  OU 

00  11 


249/1448Z-249/1445 
13A  13D 

245  47 


249/1620Z-249/1619 
7040988  9050B8D 

A071CE3  9C71 


248/1408Z-248/1406 
00  12 

2F  00 


248/1517Z-248/1406 
00  12 

2E  00 


248/1701Z-248/1406 
2000565  4010A7C 

C031A9A  47D 


248/1837Z-248/1406 
9050C8F  13091AA6 

3A204DC8  16AB 
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04744  0.2693  141.187W  2 

(  1)  2010365  200066A 

492C62CF  513D69D2 

BD  DO 

04744  0.269S  141.187W  2 

(  1)  2000360  2000362 

2E133CC1  321740C3 

BD  C7 

04744  0.288S  141.227W  2 


(  1) 


2000360 

2E133CC1 


2000362 

321740C3 


BD 

C9 

04744 

0.3Q8S  141.291W 

2 

(  1) 

OA 

OE 

100045D  1000561 

BC 

Cl 

04744 

0.334S  141.330W 

2 

(  1) 

OA 

00 

2B 

3C 

BC 

EE 

04744 

0.424S  141-486W 

2 

(  1) 

OA 

00 

2C 

37 

BB 

EE 

04744 

0.441S  141.513W 

2 

(  1) 

OA 

00 

2F 

00 

00 

00 

04745 

4.429N  142.307W 

2 

(  1) 

OA 

oc 

OA 

25 

BA 

D8 

04745 

4.429N  142.307W 

2 

(  1) 

OA 

OA 

OA 

25 

BA 

E4 

248/2329Z-248/2328 
F0A149E  1F1325B4 

673385DA  26BA 


249/0114Z-248/2328 
A030E92  1507 1EA9 

411E55CC  19AF 


249/0116Z-249/0109 


A030E92 

411E55CC 


15071EA9 

19AF 


249/0400Z-249/0356 
24  241 

1000563  49 


249/0539Z-249/0535 
00  11 

30  00 


249/1218Z-249/1213 
00  OD 

31  00 


249/1358Z-249/1353 
OA  OD 

00  00 


248/1407Z-248/1404 
03  OA 

25  00 


248/1517Z-248/1404 
03  OA 

24  00 
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04745  4.446N  142.328W  2 

(  1)  00  31 

802108C  A031394 

BA  D6 


248/1659Z-248/1655 
1000250  300066F 

C04169A  373 


04745  4.461N  142.339W  2 

(  1)  1000254  100025D 

340E44C5  3F1352CB 

BB  C4 


248/1838Z-248/1835 
7010988  15051BA8 

4F1A65D2  14A8 


04745  4.495N  142.385W  2 

(  1)  100065B  2010563 

482865CF  5B3380D6 

BE  CC 


248/2334Z-248/2330 
A050E91  1C0F26B2 

76429DDE  1FB4 


04745  4.490N  142.383W  0 

(  1)  1000254  201025F 

3F2752CA  4E2C66D1 

BF  F9 


249/01152-248/2335 
9050B8E  190F1FAD 

643D82D9  lAAF 


04745  4.526N  142.395W  2  249/01162-249/0111 

(  1)  1000254  201025F  9050B8E  190F1FAD 


3F2752CA 


4E2C66D1 


BF 

FA 

04745 

4.557N  142.425W 

2 

(  1) 

OE 

13A 

B021798  E031B9E 

BE 

FI 

04745 

4.586N  142.431W 

2 

(  1) 

OA 

OD 

OA 

25 

BD 

El 

04745 

4.645N  142.499W 

2 

(  1) 

OA 

OB 

OA 

25 

BC 

D8 

04745 

4.656N  142.522W 

2 

(  1) 

OA 

OA 

OA 

24 

BC 

C6 

643D82D9  lAAF 


249/0401Z-249/0357 
100035F  400097B 

110423A5  585 


249/0541Z-249/0537 
02  OA 

25  00 


249/1217Z-249/1211 
02  OA 

25  00 


249/1355Z-249/1352 
01  OA 

24  00 


ARGOS  READY 
/LO 
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sI^SCOTTX  EXPERIME .  NTS  \  JGOFS  \  BUOY  \  0  4  7  4  2 A .  CAL 
^sday  July  28,  1993  07:28:02  pm 

llration  File  -  OCM  Shiplaunch  ID4742  OCR-100  S/N  007  ED-100  S/N  010 
lid  for  Data  collected  between  Aug  10/92  and  Oct/ 92 
lAbration  Date  Aug  10/92  -  Lamp  S721 

K  Lu683,  Lu670,  Lu555,  Lu510,  Lu490,  Lu443,  Lu412,  Ed490 
^equation  CalRad  =  scale  *  (counts  -  offset) 

^-3  0 
0 

n-3  0 
9e-3  0 

S-3  0 

^  0 


Page:  i 


t;\USR\SCOTT\EXPERIME.NTS\JGOFS\BUOy\INFILES.LST 
Wednesday  July  26,  1993  07:42:09  pm 

AUG31.ARG 

SEPOl.ARG 

SEP02,ARG 

SEP03.ARG 

SEP04.ARG 

SEP05.ARG 

SEP06.ARG 

SEP08.ARG 

SEP09.ARG 

SEPIO.ARG 

SEPll.ARG 

SEP14.ARG 

SEP16.ARG 

SEP18.ARG 

SEP19.ARG 

SEP21.ARG 

SEPT22.ARG 

SEP24.ARG 


Page 


\  USR\ SCOTT\ EXPERIME . NTS \ JGOF S \ BUOY \ OCM_RUN 1 
dnesday  July  28,  1993  07:28:04  pm 

C:\USR\SCOTT\EXPERIME . NTS\ JG0FS\BU0Y>0CM_PR04 
ENTER  output  file  name: 

4742_05.DAT 

ENTER  platform  calibration  file  name: 
04742A.CAL 

processing  hit  1  platform  04742 
processing  hit  2  platform  04742 
processing  hit  3  platform  04742 
processing  hit  4  platform  04742 
processing  hit  5  platform  04742 
processing  hit  6  platform  04742 
processing  hit  7  platform  04742 
processing  hit  8  platform  04742 
processing  hit  9  platform  04742 
processing  hit  10  platform  04742 
processing  hit  11  platform  04742 
processing  hit  12  platform  04742 
processing  hit  13  platform  04742 
processing  hit  14  platform  04742 
232  lines  processed 
14  data  hits  processed 


Page;  i 


W:  \ USR\ SCOTT\ EXPERIME  -  NTS  \  JGOFS \  BUO  Y '  OCM_RUN2 
Wednesday  July  28,  1993  07:28:04  pm 

I  C:  \USR\SCOTT\EXPERIME.NTS\JGOFS\BUOY>OCM_PR04 

ENTER  output  file  name: 

I4745_05.DAT 

ENTER  platform  calibration  file  name: 
04745A.CAL 

processing  hit  1  platform  04745 

■processing  hit  2  platform  04745 

processing  hit  3  platform  04745 

processing  hit  4  platform  04745 

■processing  hit  5  platform  04745 

processing  hit  6  platform  04745 

processing  hit  7  platform  04745 

■processing  hit  8  platform  04745 

processing  hit  9  platform  04745 

processing  hit  10  platform  04745 

processing  hit  11  platform  04745 

1232  lines  processed 

11  data  hits  processed 
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C:  \USR\SCOTT\EXPERIMZ.NTS\JGOFS\BUOY\4742_05  .DAT 

rdnesday  July  28,  1993  07:28:03  pm 

ID  DAY  LAT  LON  ACCOR  Lu683  Lu670 


u412  ED490  Temp  Vbat  Tick 

f742  248.5785  36.748  -121.865 

0000  0.0  14.7  15.0  15.0 

04742  248.6319  36.746  -121.864 

(.0271  9.2  14.5  15.0  33.0 

4742  248.7007  36.745  -121.863 

.2440  82.8  14.5  15.0  13,0 

(4742  248.7667  36.745  -121.863 

.4570  1172.5  33.6  15.0  7.0 

4742  248.9181  36.747  -121.867 

0.5016  151.7  14.7  15.0  30.0 

14742  248.9861  36.747  -121.866 

.3660  119.5  14.7  15.0  9.0 

04742  249.0326  36.747  -121.866 

^.2576  87.4  14.7  15.0  20.0 

■4742  249.1042  36.746  -121.867 

^.0271  13.8  14.7  15.0  1.0 
04742  249.1736  36.749  -121.864 

K.5893  1098.9  0.3  10.5  57.0 

4742  249.4326  36.747  -121.863 

0.0000  0.0  14.7  15.0  0.0 

|C4742  249.5049  36.745  -121.867 

■).0000  0.0  14.7  15.0  45.0 

T;4742  249.5708  36.750  -121.868 

^.0000  0.0  14.7  15.0  21.0 

■'4742  249.6167  36.745  -121.862 

^.0000  0.0  14.7  15.0  28.0 

04742  249.6806  36.755  -121.848 

■).1356  717.3  26.7  10.5  56.0 


Lu555  Lu510  Lu490  Lu443  L 

0  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 

350  0.0000  0.0000  0.0286  0.0268  0.0365  0.0289 

150  0.0362  0.0437  0.3432  0.4024  0.4560  0.2891 

150  0.0543  0.0547  0.4862  0.5902  4.6516  3.6862 

350  0.0724  0.0766  0.6864  0.7780  0.8574  0.5782 

150  0.0634  0.0656  0.5434  0.6171  0.6567  0.4337 

150  0.0543  0.0437  0.4004  0.4293  0.4743  0.2891 

350  0.0091  0.0109  0.0572  0.0537  0.0547  0.0289 

350  0.0000  0.0000  0.0000  0.0000  4.6516  3.6862 

350  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 

350  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 

1000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 

350  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000 

0  0.0272  0.0328  0.2002  0.2415  0.2554  0.2024 


M\USR\SCOTT\EXPERIME.NTS\JGCFS\BUOY\4745_05.dat  Page:  1 

Wednesday  July  28,  1993  07:28:03  pm 

I  10  DAY  LAT  LON  ACCUR  Lu683  Lu670  Lu555  Lu510  Lu490  Lu443  L 

u412  ED490  Temp  Vbat  Tick 

t745  248.5882  4,429  -142.307  350  0.0100  0.0000  0.0000  0.0000  0.0000  0.0000 

0000  0.0  24.8  15.0  24.0 

745  248.6368  4.429  -142.307  350  0.0100  0.0000  0.0000  0.0000  0.0000  0.0000 

0.0000  0.0  24.8  15.0  36.0 


■4745 

1.1503 

248.7076  4.446 

14.9  24.8  15.0 

-142.328 

22.0 

350 

0.0100 

0.0000 

0.0269 

0.0771 

0.1439 

0.1385 

04745 

^.9896 

248.7764  4.461 

99.6  24.9  15.0 

-142.339 

4.0 

350 

0.0199 

0.0107 

0.1882 

0.5398 

0.9350 

0.8725 

■4745 

".4781 

248.9819  4.495 

154.4  25.4  15.0 

-142.385 

12.0 

350 

0.0199 

0.0214 

0.2689 

0.7197 

1.2947 

1.2603 

04745 
■  .2526 

249.0521  4.490 

129.5  25.6  15.0 

-142.383 

57.0 

0 

0.0199 

0.0214 

0.2420 

0.6426 

1.1328 

1.0803 

■4745 

0.2255 

249.0528  4.526 

89.6  -2.1  10.5 

-142.395 

18.0 

350 

0.0199 

0.0214 

0.2420 

0.6426 

0.3237 

0.2493 

■4745 
■  .2129 

249.1674  4.557 

24.9  25.4  15.0 

-142.425 

49.0 

350 

0.0100 

0.0000 

0.0269 

0.1028 

0.1978 

0.1939 

04745 

^.0000 

249.2368  4.586 

0.0  25,2  15.0 

-142.431 

33.0 

350 

0.0100 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

■4745 

%.0000 

249.5118  4.645 

0.0  25,1  15.0 

-142.499 

24.0 

350 

0.0100 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

04745 

^.0000 

249.5799  4.656 

0.0  25.1  15.0 

-142.522 

6.0 

350 

0.0100 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

I 

I 

I 

I 


